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ABSTRACT 


4 

Predictive  capability  for  the  intensity  of  severe  vortical  storms  of  the  tropo¬ 
sphere  entails  anticipation  of  which  systems  will  undergo  transition  from  one¬ 
cell  to  two-cell  structure,  and  which  will  not.  Roughly  half  of  all  tropical 
storms,  and  nearly  half  of  all  mesocyclones,  undergo  this  transition,  in  which 
the  low-level  pressure  differential  between  periphery  and  center  may  increase 
from  0(1%)  to  0(10%)  of  atmospheric  pressure,  and  the  peak  swirl  speed  may 
increase  from  0(50  m/s)  to  0(100  m/s).  The  transition  entails  insertion  of 
relatively  dry,  nonrotating,  compressionally  heated,  tropopause-level  air  at  the 
axis  of  the  vortex,  such  that  slowly  recirculating  air  (an  eye)  displaces  to  an 
annulus  (the  eyewall)  the  rapidly  swirling  cloudy  updraft  arising  from  separa¬ 
tion  of  the  low-level  inflow  layer.  This  investigation  inquires  whether  the 
vertical  stratification  of  the  circulation  at  the  periphery  of  the  vortex 
(which  varies  only  slowly  in  time)  might  not  provide  a  relatively  accessible 
observable  that  suggests  whether  two-cell  structure  can  be  sustained  by  the 
vortex.^.  To  this  end,  an  Integral-type  (transversely  averaged)  numerical  solu¬ 
tion  Is  obtained  to  an  Invlscld,  thin-layer  formulation  for  the  position, 
thickness,  wind-speed  components,  and  thermodynamic  state  of  the  separated 
inflow  layer.  The  free  layer  for  the  tornadic  mesocyclone,  of  relatively 
limited  radial  extent.  Is  described  conveniently  in  two  segments:  (1)  a  low- 
level  portion,  termed  the  turnaround,  in  which  (with  increasing  altitude)  the 
separated  layer  undergoes  an  Incursion  toward  the  axis  of  rotation  and  then  an 
excursion  away  from  the  axis,  such  that  the  layer  is  inferred  to  become  un¬ 
stable  and  undergo  a  breakdown  when  the  fluid  recovers  a  radial  distance  com¬ 
parable  to  that  at  separation;  and  (2)  a  vertically  extensive  portion  (extend¬ 
ing  to  near-tropopause  height)  in  which  the  layer  remains  of  fairly  constant 


thickness  but  is  situated  modestly  further  from  the  axis  with  increasing  alti¬ 
tude.  Peripheral  circulation  constant  with  altitude  characterizes  solutions 
compatible  with  two-cell  structure.  For  the  relatively  broad  horizontal  ex¬ 
panse  of  a  tropical  cyclone,  the  turnaround  alone  describes  virtually  the  entire 
ascent  through  the  troposphere  of  the  separated  inflow  layer;  furthermore,  the 
separated  layer  thickens  until  it  is  no  longer  treated  well  by  the  thin-strip- 
type  formulation  undertaken  here.  Nevertheless,  the  theory  does  give  insight 
into  what  balance  of  forces  describes  conservation  of  momentum  in  a  typhoon 
eyewall . 


viii 


1.  Introduction 

Progress  In  forecasting  the  Intensity  of  already  organized  vortical  storms 
requires  Identifying  observable(s)  that  foretell  whether  or  not  evolution  from 
one-cell  to  two-cell  structure  will  occur.  Very  roughly,  about  half  of  all 
tropical  storms  (Fendell  1974),  and  about  one-quarter  to  one-half  of  meso- 
cyclones  (Brooks  1949;  Houze,  Jr.  and  Hobbs  1982),  undergo  transition  from  a 
Rankine-type  radial  profile  for  swirling  (with  low-level  inflow,  central  up¬ 
draft,  and  upper-level  outflow).  Specifically,  with  the  insertion  of  a  central 
column  of  relatively  dry,  clear,  nonrotating,  slowly  recirculating,  originally 
tropopause-level  air,  the  region  of  intensely  swirling,  cloudy  updraft  becomes 
an  annulus  displaced  from  the  axis  of  symmetry.  The  full  insertion  of  an  "eye" 
within  an  "eyewall"  Implies  that  ground-level  pressure  deficits  from  the  center 
to  the  ambient  edge  of  the  vortex,  on  the  order  of  one-or-two  percent  of  atmos¬ 
pheric  pressure  (and  peak  swirl  of  not  much  over  50  m/s,  at  most)  prior  to 
insertion  of  an  eye,  increase  to  deficits  on  the  order  of  five-to-ten  percent 
of  atmospheric  pressure  (and  peak  swirl  modestly  in  excess  of  100  m/s)  (Miller 
1958).  Is  there  an  accessible  observable  that  would  permit  anticipating  whether 
a  tropical  storm  will  become  a  supertyphoon  and/or  whether  a  tornado  cyclone 
will  spawn  one  or  more  virulent  tornadoes  (M  h  duration,  150  km  path,  100  m/s 
peak  swirl)?  Of  course,  the  intensity  of  a  tropical  storm  is  not  a  simple  mono¬ 
tonic  Intensification,  followed  by  a  simple  monotonic  decay;  rather,  cyclic 
strengthening-weakening  suggests  several  partial  insertions  and  retractions 
of  an  "eye"  (Carrier  1971b);  damage  along  the  path  of  a  tornado  [and  peak  swirl 
speed  inferred  from  funnel-cloud  length  (Dergarabedian  and  Fendell  1971a, 

1971b)]  suggests  similar  behavior. 
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Accessibility  suggests  seeking  the  observable  at  the  periphery  of  the  vor¬ 
tex;  e.g.,  the  upper-tropospheric  cloud  cover  can  frustrate  passive  sensors  of 
a  meteorological  satellite.  It  is  well  known  that  a  midtropospheric  minimum  in 
the  total  static  enthalpy  (the  sum  of  static  enthalpy,  latent-heat  equivalent 
of  water-vapor  content,  and  gravitational  potential  energy)  portends  severe 
weather  (Darkow  1967);  however,  the  nature,  time  scale,  and  degree  of  severity 
usually  cannot  be  anticipated  simply  on  the  basis  of  a  local,  convectively 
unstable  stratification.  Thus,  one  turns  to  the  wind  magnitude  and  direction: 
specifically,  js.  the  variation  with  al titude  of  the  circulation  at  the  edge  of 
a^  moderate  one-cell  vortex  indicative  of  the  capacity  of  the  system  for  inten¬ 
sification  to  two-cell  structure?  It  is  anticipated  that  the  vertical  profile 
of  the  circulation  at  the  edge  of  the  vortex  (1)  varies  but  modestly  over  sys¬ 
tem  lifespan  and  (2)  varies  but  modestly  over  radial  position  (until  one 
reaches  the  central  regions,  or  "core",  of  the  vortex).  Thus,  the  stability 
of  the  vortex,  if  it  were  to  undergo  transition  to  two-cell  structure,  plausibly 
may  be  Indicated  by  one  sampling  of  variation  of  peripheral  circulation  with 
altitude. 

Application  of  linear  stability  theory  to  a  nascent  weak  disturbance  to 
identify  criteria  for  intensification  is  difficult  because  many  processes  main¬ 
taining  the  ambient  today.must  be  described  by  parameterizations  of  uncertain  valid¬ 
ity;  also,  early  entrance  of  nonlinear  effects  frustrates  identification  of  the 
state  to  which  the  system  evolves.  Applying  linear  stability  to  an  already  well- 
defined  vortical  storm  would  seem  more  relevant  for  Identifying  conditions  for 
evolution  to  two-cell  structure;  however,  the  same  difficulties  with  parameteri¬ 
zations  (can  sub-vortex-scale  phenomena  such  as  •  imulus  convection,  turbulent 


diffusion  and,  possibly  relevantly,  radiative  transfer  be  parameterized  in  terms  of 
vortex-scale  variables  and,  if  so,  what  are  the  correct  parameterizations?)  and 
with  smal  1  -perturbation  theory  hold- -and  adequate,  yet  simple,  description  of  a  one- 
cell  configuration,  for  facile  perturbation,  is  not  a  trivial  matter.  Instead 
of  a  temporal  development,  a  steady  analysis  is  undertaken  that  inquires  whether 
"eye-within-an-eyewall"  configuration  is  consistent  with  only  certain  ambient 
stratifications  of  circulation.  If  it  turns  out  that  two-cell  structure  is 
compatible  with  virtually  any  variation  of  ambient  circulation  with  altitude, 
then  no  discriminant  has  been  identified.  Also,  it  must  turn  out  (1)  that  the 
analysis  is  insensitive  to  parameterization  details,  or  else  it  suffers  from 
limitations  already  noted,  and  (2)  that  results  are  insensitive  to  geometric 
details,  or  else  the  approximate  analysis  to  be  undertaken  is  too  crude. 

Finally,  it  is  remarked  that  the  translational  path  of  the  vortex  is  an 
entirely  separate  matter  that  can  be  approached  by  statistical  data  from  histori¬ 
cal  record,  or  by  embedding  the  vortex  in  larger-scale  atmospheric  dynamics. 
Attention  here  is  confined  to  anticipating  intensity. 

2.  Structure  of  a  Severe  Atmospheric  Vortex 
2.1  Four-Part  Structure  of  a  Two-Cell  Vortex 

Since  the  possible  transition  from  one-cell  to  two-cell  structure  primarily 
concerns  alteration  in  the  more  rapidly  swirling  "core"  of  the  vortex,  a  time- 
averaged  axlsymmetric  analysis  Is  undertaken.  It  seems  premature  to  account  for 
gustiness  and  asymmetries  [e.g.,  manifestation  of  large-scale  coherent  struc¬ 
ture  (Roshko  1976)  such  as  suction  vortices  (Fujita  et  al  1976)]  until  a  more 


i  i. 


elementary  description  Is  established.  It  is  taken  that  time-dependent  investi¬ 
gation  of  intensification  and  of  decay  would  yield  results  consistent  with  the 
structure  proposed  for  a  postulated  quasisteady  mature  stage. 

A  closed  model  is  adopted,  such  that  diffusive  (but  not  convective)  trans¬ 
port  across  domain  boundaries  may  occur.  These  boundaries  are  a  lateral  edge  at 
which  the  swirl  has  diminished  to  virtually  ambient  value;  an  isobaric,  isothermal 
slippery  1  id  set  at  the  altitude  to  which  ground-  (or  sea-)  level  air  would  rise, 
in  view  of  its  total  static  enthalpy  (in  the  convectively  unstable  stratification), 
upon  moist-adiabatic  ascent^;  and  the  ground  plane  or  sea  surface,  approximated 
as  an  Isothermal  boundary  enforcing  a  virtual  no-slip  constraint  on  the  velocity 
field.  In  a  numerical  approach,  a  uniformly  valid  boundary-value  problem  would 
be  posed  and  solved  for  the  postulated  quasisteady  mature  stage  of  a  severe 


*In  the  locus  of  thermodynamic  states  referred  to  as  the  moist  adiabat,  air 
expands  and  cools  as  it  rises,  until  it  becomes  saturated.  Thenceforth 
during  ascent,  sufficient  water  vapor  condenses  to  maintain  saturation  at  the 


issa 


locally  prevailing  temperature.  The  air  retains  the  heat  of  phase  transition, 
and  the  condensed  water  substance  either  falls  out  or  is  conveyed  along,  but 
reversal  of  the  process  (evaporative  cooling)  is  precluded.  Ascent  is  idealized 
as  too  rapid  to  permit  dilution  of  the  updraft  with  drier  and/or  cooler  air.  The 
level  of  neutral  stability  entails  recovery  of  ground-level  (or  sea-level)  total 
static  enthalpy;  identification  of  the  system  "lid"  at  such  an  altitude  deviates 
from  conventional  definition  of  the  tropopause.  In  fact,  of  course,  in  idealiz¬ 
ing  the  eyewall  as  characterized  by  a  moist  adiabat,  one  is  ignoring  structure 
within  the  eyewall,  and  Is  omitting  dilution  of  eyewall  air  with  air  from  the 
eye  and/or  with  air  from  the  potential  vortex  (Appendix  A). 
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atmospheric  vortex.  In  the  subdivisional  approximate  analyses  to  be  undertaken 
here,  only  the  locally  dominant  processes  are  retained,  and  enforcing  appro¬ 
priate  continuity  of  fields  and  fluxes  at  subdivisional  interfaces  yields  a 
global  solution;  self consistency  is  established  by  a  posteriori  substitution 
of  the  solution  to  confirm  that  the  locally  dominant  processes  have  been  anti¬ 
cipated  correctly.  It  is  recognized  that  some  relaxation  of  the  closed  nature 
of  the  model  may  be  required  (e.g.,  possible  downflux  of  stratospheric  air  into 
an  eye,  upper- tropospheric  outflow  of  eyewall  air,  etc.};  nevertheless,  the 
closed,  subdivisional  model  is  a  convenient  vehicle  for  presentation. 

A  previously  adopted  three-  or  four-part  structure  for  a  mature,  quasisteady 
axisymmetric  severe  atmospheric  vortex  (Carrier  1970;  Carrier  et  al .  1971)  is 
briefly  reviewed  to  permit  the  presentation  to  be  sel fcontained  (figure  1). 

The  bulk  of  the  vortex  consists  of  rapidly  swirling  air  (region  I),  the  radial 
variation  of  the  azimuthal  velocity  component  being  approximated  adequately 
as  that  of  a  potential  vortex;  the  ultimate  source  of  the  swirl  is  the  rota¬ 
tion  of  the  earth.  Whereas  there  Is  appreciable  radial  influx  in  the  bulk  of 
the  vortex  during  spin-up,  in  the  quasisteady  mature  stage  of  interest  here, 
the  radial  Influx  need  preclude  only  outward  movement  of  the  central -ascent 
region,  and  since  only  a  very  small  influx  Is  required  for  that,  the  radial 
Influx  seems  negligible.  Thus,  after  spi/i-up,  there  is  a  cyclostrophic  bal¬ 
ance  [holding  radially  between  the  central  updraft  and  peripheral  edge,  hold¬ 
ing  axially  between  the  surface-layer  inflow  and  upper-troposphere  outflow  (see 
below)],  such  that  the  radial  pressure  gradient  balances  the  centrifugal  force; 
the  more  meticulous  (so-called)  gradient-wind  balance  does  not  seem  required 


for  present  purposes.  The  angular  momentum  at  any  radial  position  in  the  bulk 
of  the  vortex  is  that  holding  at  the  periphery  at  the  same  height.  Whereas 
there  may  be  local  transient  updrafts  and  downdrafts,  on  average  there 
is  a  slow  downdraft  of  air  from  the  bulk  vortex  to  the  surface  inflow  layer; 
while  the  local  downflux  may  be  modest  in  magnitude,  the  total  downflux,  upon 
summing  over  area,  suffices  to  furnish  the  "throughput"  of  the  vortex  (i.e., 
the  secondary  flow,  associated  with  low-level  influx,  central  updraft,  and 
upper-level  outflow).  The  thermodynamic  fields  throughout  the  "throughput- 
supply"  region  (i.e.,  the  bulk-vortex  region)  are  those  prescribed  to  hold  at 
the  peripheral  edge;  the  major  exception  is  the  Isobars,  which  are  horizontal 
at  the  periphery  (where  hydrostatics  holds),  but  which  dip  groundward  (or  sea¬ 
ward)  with  enhanced  swirling  (I.e.,  with  decreasing  radial  distance  from  the 
axis  of  symmetry).  While  the  thermodynamic  stratification  of  the  bulk-vortex 
region  is  Indicative  of  the  ambient  in  which  the  vortex  was  formed,  sea-to-alr 
transfer  maintains  the  stratification  for  the  tropical  cyclone  over  the  trades 
in  a  manner  precluded  for  a  continental  mesocyclone.  In  any  case,  the  through¬ 
put  supply  Is  finite,  and  the  later  throughput  may  be  drier  and  cooler  than 
earlier  throughput;  however,  other  phenomena  may  enter  to  disrupt  the  system 
prior  to  throughput  exhaustion  (e.g.,  disruption  by  orographic  lifting). 

The  low-level  inflow  region  (II)  is  the  only  one  in  which  angular  momentum  is 
not  conserved;  it  is  dissipated  in  part,  owing  to  the  no-slip  boundary  condition 
that  engenders  a  swirling  Influx  (so-called  "tea-cup  effect").  Thus,  indeed  there 
must  be  a  modest  downdrift  on  average  into  the  surface  layer,  in  which  radial 
and  azimuthal  velocity  components  are  comparable  and  significant.  In  fact,  there  is 


an  accelerating  radial  pressure  gradient  acting  on  the  inflow  layer;  thus  it  is 
not  surprising  that,  under  the  high-speed  portion  of  the  vortex,  the  diffusive- 
like  processes  are  confined  to  a  thin  sublayer  contiguous  to  the  ground  or  sea-- 
a  sublayer  that  becomes  ever  thinner  as  the  axis  of  rotation  is  approached,  even 
though  the  total  inflow-layer  thickness  grows  (albeit  modestly)  with  decreasing 
radial  distance  from  the  axis  (Burggraf  et  al  1971;  Carrier  1971a;  Carrier  and 
Fendell  1978).  Specifically,  over  the  bulk  of  the  inflow  layer,  inviscid  pro¬ 
cesses  are  controlling;  at  fixed  radial  distance  from  the  axis,  as  one  descends 
vertically  downward,  the  virtual  pure  swirling  of  the  bulk  vortex  gradually 
becomes  virtual  pure  radial  influx  of  the  same  magnitude  (appropriate  for  the 
radial  distance  from  the  axis).  Only  across  the  sublayer  does  the  velocity  fall 
abruptly  to  zero.  This  established  behavior  is  belabored  because  (1)  it  justi¬ 
fies  an  inviscid  treatment  of  most  of  the  Inflow  layer  if  it  separates  to  form 
an  eyewall;  (2)  it  emphasizes  that  the  portion  of  the  eyewall  contiguous  to  the 
eye  may  shear  the  eye  into  recirculation,  but  can  impart  little  rotation  to  the 
eye  fluid;  (3)  it  emphasizes  that  the  portion  of  the  eyewall  contiguous  to  the 
potential  vortex  but  modestly  shears  the  potential  vortex;  (4)  it  emphasizes  the 
nature  of  the  (sublayer)  fluid  that  lies  below  the  eye,  _1f  the  eye  does  not 
extend  from  the  tropopause  down  the  full  extent  of  the  troposphere  to  the  sea 
surface  (or  ground);  and  (5)  It  characterizes  the  velocity-component  profiles 
holding  In  the  inflow  layer  at  separation. 

Whereas  the  quantitative  description  of  the  bulk  vortex  and  surface- Inflow 
layer  seems  already  adequately  in  hand  for  present  purposes,  the  description  of 
the  turnaround  [whereby  the  separating  inflow  layer  turns  vertical  (lowest 


portion  of  III)]  and  the  eyewall  [whereby  the  vertical  inflow  layer  rises  the 
extent  of  the  troposphere  (middle  and  upper  portion  of  III)]  has  been  but  very 
roughly  examined,  and  Is  the  subject  of  closer  scrutiny  here.  It  may  be  noted 
that  for  a  one-cell  vortex,  the  surface-inflow  layer  continues  into  the  axis 
(r*=0)  without  separation,  and  then  turns  upward  by  necessity.  The  emphasis  on 
inflow-layer  separation  into  a  free  shear  layer  (the  position  of  whose  inner  boun¬ 
dary/interface  with  the  eye,  and  the  position  of  whose  outer  boundary/ interface 
with  the  bulk  vortex,  are  set  by  requiring  continuity  of  pressure  and  are  ascer¬ 
tained  in  the  course  of  solution)  anticipates  that  two-cell  structure  is  to  be 
scrutinized  below.  It  may  also  be  noted  that  a  moist-adiabatic  locus  of  thermo¬ 
dynamic  states  holds  in  the  separated  layer,  that  angular  momentum  Is  conserved 
along  its  streamlines,  that  the  total  head  along  its  streamlines  is  constant  (aside 
from  the  gravitational -potential  term),  that  a  cyclostrophic-like  force  balance 
holds  In  the  layer,  and  that  the  separating-inflow-layer  thickness  and  the  turn- 
around-and-eyewall  thickness  are  probably  roughly  comparable  (at  least  up  to 
midtroposphere).  Thus,  the  turnaround  and  eyewall  Involve  a  relatively  thin 
layer  across  which  the  transverse  pressure  gradient  is  significant;  frequently 
In  highly  subsonic  flow,  the  transverse  pressure  gradient  is  small  across  thin 
layers,  so  the  present  context  is  not  entirely  conventional,  and  fortunately 
Integral -type  (transversely  averaged)  treatment  is  adequate  for  present  pur¬ 
poses.  Whereas  the  upper-tropospheric  turning  of  the  eyewall  into  an  outflow 
layer  Is  to  be  the  subject  of  closer  scrutiny  in  future  work,  here  attention  is 
concentrated  on  the  lower  troposphere  and  midtroposphere. 

The  locus  of  thermodynamic  states  that  characterize  the  inflow  layer  and 
turnaround/eyewall  Is  given  by  the  moist  adiabat  based  on  ground-level  (or  sea- 


level)  ambient;  the  ambient  altitude  is  assigned  to  the  level  of  neutral  stability. 
The  eye  (IV)  is  idealized  as  completely  dry,  so  its  locus  of  thermo¬ 
dynamic  states  is  based  on  a  dry  adiabat,  referenced  to  the  neutral-stability-lid 
altitude  and  state.  Mass  exchange,  from  the  eye  to  the  eyewall  at  the  lower  alti¬ 
tude  and  from  the  eyewall  to  the  eye  at  greater  altitude,  is  to  be  included  in 
future  investigations;  incidentally,  the  net  eye/eyewall  exchange  need  not  be  in 
perfect  balance,  though  imbalance  implies  net  flux  across  the  ceiling  lid.  It  is 
in  the  emphasis  on  the  significance  of  the  eye  for  the  intense  vortex  that  the 
structure  presented  here  is  distinguished  from  other  treatments  (e.g.,  Barcilon 
1967). 

2.2  Implications  of  Lateral  Scale  on  Vortex-Core  Analysis 

The  severe  atmospheric  vortex  has  been  characterized  as  a  once-through 

2 

cycle,  in  which  a  throughput  supply  is  exhausted  as  air  is  processed  ;  the 
processed  air  is  too  dry  for  effective  recycling,  and  recycling  on  the  scale  of 
vortex  lifetime  might  entail  compressional  heating  (during  descent  at  the  periph¬ 
ery)  incompatible  with  a  central -pressure  deficit  at  ground  level  (or  sea  level). 

As  already  noted,  specification  of  the  ambient  stratification  (holding  at 
r*=r0*)  implies  the  amount  of  lateral ,  low-altitude  reduction  of  pressure  at  the 
axis  of  symmetry  from  the  value  of  pressure  holding  at  the  periphery,  on  the 
basis  of  moist  adiabatic  ascent  (one-cell  vortex),  and  the  amount  of  further 
pressure  reduction  on  the  basis  of  dry-adiabatic  compression  of  air  from  the 
neutral -stability  level  (two-cell  vortex)  (Lilly  1969;  Gray  1970;  Dergarabedian 

The  closed  model  with  a  "captured"  throughput  supply  is  more  appropriate  for 
the  large  lateral  scale  of  the  tropical  cyclone  than  for  the  limited  lateral 
scale  of  a  rotating  thunderstorm. 
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and  Fendell  1970,  1977).  The  eye  air  need  not  be  recompressed  eyewall  air,  but 
could  be  descending  stratospheric  air  (partly  or  entirely).  Contrary  to  postu¬ 
lations  of  earlier  models  (Malkus  and  Riehl  1960),  recent  satellite  photography 
(Black  1977)  confirms  sufficient  outward  slope  of  the  eye/ eyewall  interface 
(Haurwitz  1935),  so  that  the  augmented  pressure  deficit  from  ambient  (associated 
with  eye  insertion)  is  available  to  support  enhanced  eyewall  swirl  speeds.  In 
particular,  if  one  cencentrates  attention  just  above  the  surface  inflow  layer, 
ignores  lateral  density  variations  at  low  altitude  in  the  dynamics,  and  adopts  a 
radial  profile  for  the  swirl,  integration  of  the  cyclostrophic  approximation 
to  the  conservation  of  radial  momentum  yields  the  peak  swirl  for  a  thermo- 
hydrostatically  computed  low-level  pressure  deficit.  For  example,  a  Rankine- 
vortex  model,  such  that  the  pressure  deficit  is  expended  in  equal  amounts  to 
sustain  the  rigidly  rotating  core  and  the  outer  potential  vortex,  is  plausible 
for  a  one-cell  vortex;  since  the  eye  is  nonrotating  in  two-cell  structure,  the 
entire  (larger)  pressure  deficit  sustains  an  outer-potential -vortex  profile 
(only).  In  any  case,  one  obtains  a  peak  swirl  speed  holding  near  ground  level 
at  the  innermost  radius  of  the  potential  vortex,  V*.  As  already  noted,  V*  is 
about  50  m/s  or  less  for  one-cell  structure  (moist-adiabatic  ascent  near  the 
axis),  and  over  100  m/s  for  two-cell  structure  (dry-adiabatic  eye  near  the 
axis)— since  tornado-prone  and  hurricane-prone  ambients  often  yield  comparable 
low-level  central  pressure  deficits  from  ambient. 

The  other  parameter  that  characterizes  a  vortex  Is  r-j*,  the  radial  dis¬ 
tance  from  the  axis  at  which  V*  holds;  also  at  r*=r-|*,  the  inflow  layer  (II) 
separates  to  constitute  the  turnaround  (lowest  portion  of  III)  for  cases  of 


two-cell  structure  (i.e.,  cases  in  which  region  IV  exists).  For  an  intense  tor¬ 
nado,  this  distance  is  taken  to  be  0(170  m),  while  for  a  typhoon  this  distance 
is  taken  to  be  0(16  km).  From  knowledge  of  V*  and  r^*,  and  of  the  circulation 
and  thermodynamic  state  as  a  function  of  altitude  at  the  vortex  periphery  (where 

ro*>>rl*’  so  t*1at  t^ie  PeriPhery  may  be  taken  to  be  at  infinity,  without  sig¬ 
nificant  error,  for  many  purposes),  one  expects  to  be  able  to  deduce  all  other 
mean  properties  of  the  vortex. 

Clearly,  both  tornado  and  hurricane  extend  through  the  troposphere,  or 
0(14  km)  vertically;  however,  the  vast  difference  in  radial  scale  between  a 
devastating  tornado  arising  from  evolution  to  two-cell  structure  in  a  meso- 
cyclone  and  a  devastating  typhoon  from  evolution  to  two-cell  structure  in  a 
tropical  cyclone  necessitates  special  provision  for  treatment  of  the  two  vor¬ 
tices  within  one  formulation.  In  particular,  it  Is  anticipated  that  the  sepa¬ 
rating  Inflow  layer  (perhaps  50  m  thick  for  a  tornado,  and  perhaps  1  km  thick 
for  a  typhoon)  continues  radially  inward  as  it  rises  after  separation,  reaches 
a  minimum  finite  distance  from  the  axis,  and  then  turns  outward  from  the  axis. 

It  is  further  anticipated  that  the  altitude  over  which  this  inward,  outward 
(turnaround)  excursion  occurs  is  0(rj*).  For  a  tornado,  the  entire  turnaround 
involves  the  lowest  170  m  or  so  of  the  troposphere.  The  overwhelming  prepon¬ 
derance  of  the  description  of  the  fate  of  the  separated  inflow  layer  lies  sub¬ 
sequent  to  the  turnaround,  in  what  was  termed  the  eyewall.  For  a  typhoon,  the 
turnaround  might  involve  a  vertical  distance  of  14  km  or  so,  the  entire  depth 
of  the  trophsphere;  if  the  entire  fate  of  the  separated  inflow  layer  is  given 
by  the  turnaround,  then  there  is  no  sequential  eyewall  region  after  separation. 


2.3  The  Turnaround  for  the  Typhoon  and  for  the  Hurricane 

For  the  typhoon,  one  anticipates  the  need  to  retain  the  gravitational 
potential,  the  variation  of  the  bounding  eye  and  eyewall  pressures  with  alti¬ 
tude,  and  the  change  of  density  with  altitude,  in  treatment  of  the  turnaround. 
However,  for  the  tornado,  one  may  ignore  change  of  the  gravitational  poten¬ 
tial,  change  of  the  bounding  eye  and  eyewall  pressures  with  altitude,  and  the 
change  of  density  with  height;  indeed,  one  may  ignore  even  the  variation  of 
thickness  of  the  separated  layer  with  height  for  the  tornado,  since  the  thick¬ 
ness  is  so  small.  Thus,  a  simple,  semianalytic  incompressible  model  may  be 
derived  for  the  tornado,  since  its  only  purpose  is  to  furnish  conditions  hold¬ 
ing  at  the  end  of  the  inward-outward  excursion  of  the  separated  inflow  layer 
in  the  turnaround;  such  conditions  hold  at  the  start  of  the  more  vertical 
rise  of  the  eyewall.  However,  it  turns  out  that  analysis  of  the  turnaround 
for  the  tornado  may  be  superfluous  for  even  this  objective.  The  reason  is 
that  completion  of  the  inward-outward  excursion  of  the  turnaround  leads  to 
a  configuration  such  that  the  streamlines  of  the  separated  layer  are  concave 
in  the  direction  of  increasing  streamwise  velocity.  The  associated  instability 
is  taken  here  to  imply  occurrence  of  the  still  incompletely  understood  vortex- 
breakdown  phenomenon.  Thus,  no  attempt  is  made  to  infer  conditions  at  the 
outset  of  the  tornado  eyewall;  rather,  plausible  conditions  holding  after  the 
breakdown  are  estimated. 

The  density  decreases  with  altitude,  so  by  simple  continuity  considera¬ 
tions,  one  anticipates  that  the  displacement  of  the  eyewall  from  the  axis  of 
symmetry  increases  with  height;  however,  the  multiple  of  density  decrease 
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limits  the  multiple  of  increase.  While  the  direct  procedure  for  investigat¬ 
ing  the  eyewall  would  be  to  specify  the  stratification  of  ambient  circulation 
and  to  deduce  the  radial  displacement  of  the  eyewall  from  the  axis,  it  turns 
out  more  tractable  to  solve  the  inverse  problem  in  which  the  displacement  is 
specified  and  the  associated  stratification  of  circulation  is  inferred. 


3.  The  Turnaround 
3.1  The  Tornado 

An  Invlscld  incompressible  steady  axisymmetric  flow  with  swirl  could  be 
adopted  to  describe  the  lowest  portion  of  region  III;  in  general,  a  quasil inear 
second-order  elliptical  partial  differential  equation  for  the  streamfunction 
would  satisfy  nonlinear  gradient-type  constraints  on  both  the  streamsurface 
contiguous  to  the  eye  (IV)  and  also  the  streamsurface  contiguous  to  the  outer 
potential  vortex  (I)  (see,  e.g.,  Batchelor  1967,  pp.  543-546).  Such  a  chal¬ 
lenging  free-boundaries  problem  would  require  numerical  solution;  even  then,  it 
would  not  be  of  particular  value,  since  retention  of  density  variation  seems 
essential  to  any  treatment  that  aspires  to  describe  region  III  through  the  ex¬ 
tent  of  the  troposphere. 

A  transversely  averaged,  integral -type  treatment  of  the  variable-density 
separated  shear  layer  is  developed  In  succeeding  sections.  Here  a  constant- 
density  treatment  of  the  lowest  portion  of  region  III  Is  presented.  As  noted 
in  Section  2,  such  a  treatment  can  pertain  only  to  the  tornado,  in  which  the 
turnaround  entails  a  region  that  extends  from  ground  level  to  but  170  m  or  so 
(such  that  hydrostatic  effects,  and  variation  of  ambient  circulation,  are  negli¬ 
gible  over  the  vertical  extent  of  Interest).  Furthermore,  the  turnaround  for 
the  tornado  evolves  to  an  unstable  state  that  does  not  furnish  definitive  guid¬ 
ance  for  the  eyewall  calculation  (the  solution  for  region  III  over  the  rest  of 
the  troposphere).  Thus,  the  constant-density  treatment  that  immediately  follows 
serves  mainly  as  a  preparation  for  the  variable-density  treatment  required  for 
the  typhoon  turnaround  (that  occurs  over  a  height  that  extends  many  kilometers 
up  from  the  sea  surface);  the  reader  familiar  with  turnarounds  needs  but  skim 
quickly  over  the  remainder  of  Section  3.1  for  notation. 


14 


The  momentum  balance  in  an  axi symmetric,  separated  boundary  layer  without 
structure  is  sought  (  figure  2  ).  If  the  lowest  portion  of  a  fully  inserted  eye 


OH 


s 


is  isobaric  at  pressure  P*^,e,  then  the  pressure  in  the  potential  vortex  is 
(density  p*  const.) 


P*(r*,z*)  =  p*  +  4-  p*  V*s  1  -  -pi 


where  super  asterisk  denotes  a  dimensional  quantity;  at  the  separation  radius, 
p*(r1*,z*)  =  p|  ,  in  the  absence  of  hydrostatic  effects.  For  r*  <  r-|*,  the 
radial  pressure  gradient  adds  to  the  centrifugal  force  associated  with  swirl¬ 
ing  v*,  such  that  the  sum  of  two  forces  is  balanced  witn  a  centrifugal  force 
associated  with  the  (radial  and  axial)  motion  q*;  thus  q*  is  the  velocity  com¬ 
ponent  In  a  plane  containing  the  axis  of  symmetry,  and  is  perpendicular  to  v*. 

Accordingly,  if  R*(z*)  denotes  the  radial  distance  from  the  axis  to  the 
outer  surface  of  the  eyewall,  at  any  point  [R*(z*),z*],  the  force  balance 
perpendicular  to  the  thin  sheet  gives 

.  _JL  P*  =  q*2  R*" _ _ 

’  3h*  'p  [1  +  (r*')2]3/2  "  R*  [1  +  (R*')2]Vz 

where  h*  represents  a  coordinate  running  across  the  sheet,  prime  denotes  ordi¬ 
nate  derivative,  and  the  two  principal  radii  of  curvature  have  been  introduced. 


(3.1.1) 


(3.1.2) 


If  r*  =  V*  =  R*v*, 


v*2dh*  =  r*2 


[R*(h*)]: 


rV  V*2  A*  A* 

-  if  = 


R*2  [R*( h*)]2 


(3.1.3) 


’  r  t2 

q*2  dh*  ,  v*2  B*  =»  B*  =  dh*  < 


(3.1.4) 


From  these  definitions  (A*,B*  consts.),  and  from  (3.1.1)  applied  at  r* 
(3.1 .2)  integrates  to 


r^2)  A*  V*2  r-|*2 

R*2  R*3[l  +  (R*'2)]1 


B*  V*2  R*" 
[1  +  (R*')2]' 


(3.1.5) 


The  following  nondimensional ization  is  introduced: 


d  =  z  s 

V  >-1* 


z*  .  _  _h^  _  2A^  _  2B* 

r.|*  »  n  ~  r^  ’  a  ~  r-|*  *  ^  ~  ri* 


(3.1 .6) 


If  A*  ■  h*,  B*  =  (h*/2),  seemingly  reasonable  values,  then  (a/2)  =  3  =  h,  where 
h  is  0(0.2).  In  any  case,  under  (3.1.6),  (3.1.5)  becomes  (with  a, 3  specified) 


1  1  .  _ a _ BR" 

R2  R3(l  +  R'2)1/z  (1  +  R'2)3A 


r»  .  -g-p  -l 

3  R3  ft  R2  j  j 


3  R 


(3.1.7) 


For  this  translationally  Invariant  (in  z)  equation,  which  yields  R  periodic  in 
z  (see  below),  solution  is  sought  for  positive  and  negative  z,  where  the  boun¬ 
dary  conditions  are  taken  to  be 


R(0)  =  1, 


R'(0)  =  R  '  ,  given  const  (taken  positive), 


(3.1.8) 


(3.1.9) 


The  boundary  conditions  preclude  odd  or  even  solution  for  R(z).  Sought  are 
R+,  the  largest  value  of  R,  which  occurs  at  R,s  0,  and  R",  the  smallest  value 
of  R,  which  also  occurs  at  R'=  0.  Since  R"<  1,  v*(R*~)/v*(r^*)  =  [v*(R*")/V*]  >  1; 
the  magnitude  of  this  "overshoot"  of  swirl  Is  sought  as  a  function  of  (plausi¬ 
ble  values  for)  the  parameters  a, 3. 


'v  .y.v.w  •.  vv  -'v*!.  ^  .-.v.v.a' 


**  j 


While  numerical  treatment  of  (3.1.6)  -  (3.1.9)  ultimately  is  used,  some 

preliminary  phase-plane  analysis  is  useful  (figure  3): 

dR  _  n/n\ ^  d 2 R  _  _dP_  _  dR  _d£_  _  D  _dP_ 

dz  'K'  dz1  dz  "  dz  dR  "  v  dR  ’ 


hence 


0PP1  =  ct  _ 

(1  +  P2)%  R3(l  +  P2)V2  02 


+  —  -  1  . 
R2 


(3.1 .10) 


(3.1 .11) 


Although  the  following  development  is  not  pursued  here  to  the  extent  of  obtain¬ 
ing  results,  it  may  be  worth  noting  that  if 


then  (3.1.11)  becomes 


H(R)  =  (1  +  P2)‘1/z, 


-BH*  =  —  +  —  -  1  . 

R3  R2 


t  =  R-2  , 


2B  .  aH  .  T-V2  -  , 


H  •  -  (&]  !  (l  ♦  -f-J  exp(^-)  erfc  (•§■)  *  ♦  E  «p(^ 

where  E  is  a  const,  of  integration.  Since  P  =  (dR/dz), 

AD  (  -1V2 


•  {l  ♦  [H(t)]-2}  !  *  ^ 

41- - 

f  .I/2 


{l  +  [HIT)]'2- 


where  H(x)  is  given  above. 


■  *  ,**  •  * » *  ,  *  ,*»  *»  **  *»  *,  *,  *»  \  \  *•*  1 
'  *  *  *  '  1*  *  *  *  *  *  *  •  ’  O  •  *  «  *  K  *  «_«  I  1  .  1*.  t*.  4*.  A  **_  •*«.  **-  .  r 
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The  slope  P'  is  infinite  at  P  =  0,  and  is  zero  for  R  satisfying 


_ a 

R3(l  +  Pz) 


77 


(3.1.12) 


The  intersection  of  the  curve  of  (3.1.12)  with  the  ray  P  =  0  is  given  by 


_°L_  +  J_  .  i 

r*3  r77 


(3.1 .13) 


where  R*  =  1  for  a  =  0  and  R*  >  1  for  a  >  0.  For  1  >  >  a  >  0, 

R*  =  1  +  (a/2)  -  (3a2/8)  +  .  .  .  . 

For  the  special  case  of  no  swirl  a  =  0,  multiplication  of  (3.1.11)  by 
R'  yields  the  integral 

P  .  (R  -  I)2  +  f 
(1  +  R'2)1^2  R  m 


(3.1.14) 


where  the  constant  of  integration  has  been  written  as  (P/m),  with 

■  ■=  (l+R0’V/i 


(3.1.15) 


for  consistency  with  the  boundary  conditions  (3.1.8)  and  (3.1.9).  Since,  from 

V  V 

(3.1.14),  R'  is  maximum  at  R  =  1 ,  (1  +  RQ'2)  2  is  the  maximum  value  of  (1  +  R'2)  , 

whence  the  symbol  m.  At  Rf  =  0, 


R1  =  1  +  -f- 


1  - 


- 

Vz 

- 

1  ] 

2  +  -f- 

m  J 

2 

(3.1 .16) 


For  R0 '  -*■  0  so  m  -*■  1 ,  R*  merge  to  unity;  this  case  involves  a  vertically  separat¬ 
ing  surface  inflow  layer,  and  hence  no  overshoot.  For  R^  -*■«  so  m  -►  <*>,  R“  remain 


bounded : 


R*  =  1  +  -|-  ± 


6 

2 


(3.1.17) 


this  case  involves  effectively  horizontal  inflow  of  the  separating  surface  layer 
and  leads  to  the  minimum  value  of  R~  for  fixed  8.  For  3  =  0.1,  R+  =  1.37  and 
R~  =  0.73-,  for  3  =  0.2,  R+  =  1.56  and  R”  =  0.64.  Since  6  =  0(0.2)  is  a  plausible 
value,  these  values  suggest  what  proves  to  be  a  general  trend:  R"  decreases  as 

3  increases  and  as  R  '  increases. 

0 

Finite  values  for  a  indicate  finite  swirl;  increasing  a  yields  larger  values 
of  R~  and  hence  smaller  overshoot  of  swirl.  From  results  of  numerical  integra¬ 
tion  given  in  table  1  and  supplemented  by  figures  4  and  5,  for  finite  a,  R1 
approach  finite  values  as  Rq'  -*■  ®  for  fixed  3.  For  the  plausible  values  a  =  0.1, 

3  ■  0.2,  for  Rq '  =  2,  5,  10,  the  corresponding  values  of  R"  are  0.748,  0.695,  0.677. 
Hence,  swirl  overshoots  in  the  range  of  about  30%  seem  plausible,  but  not  much 
more.  Estimates  (Lewellen  1977,  p.  109)  that  the  swirl  in  the  turnaround  can 
exceed  swirl  at  R  =  1 ,  z  =  h  by  100%  are  excessive  according  to  this  analysis. 

As  suggested  earlier,  the  periodic  nature  of  the  mathematical  solution  is 
not  observed  physically  because  of  known  instability  of  flows  under  centrifugal 
forces  acting  on  concave  streamlines.  Soon  after  one  incursion  and  return  of 
the  separated  shear  layer,  "vortex  breakdown"  (Hall  1972;  Leibovich  1978)  enters, 
and  no  attempt  is  undertaken  here  to  describe  details  of  this  process.  What 
emerges  is  taken  to  be  a  vertical  free  layer,  the  evolution  of  which  with  height 
entails  a  more  Inclusive  model  (holding  through  a  considerable  depth  of  the  tro¬ 
posphere)  attempted  in  Section  4  below. 


3.2  The  Typhoon 


As  already  noted,  if  the  vertical  extent  of  the  turnaround  for  the  two¬ 
cell  tropical -storm  case  is  0{r-j*)  where  r^*  is  the  radius  of  peak  swirl  in 
the  potential  vortex  (region  I),  then  accounting  for  axial  stratification  in 
the  eye  and  in  the  ambient,  and  accounting  for  the  gravitational  potential, 
are  appropriate.  The  generalization  of  the  presentation  of  Section  3.1  to  be 
now  undertaken  is  pertinent  not  only  to  the  tropical -cyclone  turnaround;  while 
the  initial  conditions  holding  at  the  starting  value  of  the  axial  coordinate  z* 
may  be  altered,  and  while  the  roles  of  some  prescribed  and  deduced  variables 
may  be  interchanged  for  more  efficient  solution  (inverse-problem  methodology), 
still  the  ordinary  differential  equations  to  be  derived  are  as  pertinent  to  the 
eyewall  (mid  and  upper  portion  of  region  III)  as  to  the  "compressible  turnaround" 
(lower  portion  of  region  III). 

First,  some  properties  of  the  pressure  field  in  region  I  are  required;  this 
field  is  denoted  p-j*(r*,z*),  where  p^*(rQ*,z*)  =  Pi^U*).  given,  where  "amb" 
abbreviates  ambient.  Also  p1*[R*(z*) ,z*3  =  p*e[R*(z*),z*],  where  "ae"  denotes 
"^mbient-side  edge  of  the  separated  layer."  If,  for  brevity,  the  subscript 
unity  is  discarded  temporarily,  then  in  region  I,  with  subscripts  r*  and  z* 
denoting  partial  differentiation. 


Pjr*  i  p*  v*z/r*  > 


-p%  =  p*  g* 


(3.2.1) 

(3.2.2) 


hence. 


g*  pU  s 


(3.2.3) 


If  one  seeks  solution  in  the  form 


p*(n*)  *  P*[m*(z*)  -  s*(r*)J 


(3.2.4) 
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P**(n*)  |g*  s**(r*)  r*3  -  r*2(z*)  m**(z*)j  = 


(3.2.5) 


i.e.,  without  loss  of  generality, 


s*(r*)  =  -  1/(2  g*  r*2)  , 

rZ* 

m*(z*)  =  r*“2(z1)  dz1  . 

i 


(3.2.6) 


(3.2.7) 


For  the  special  case  of  a  vortex  that  is  invariant  with  altitude, 

T*  *  rQ*,  const.,  m*(z*)  =  z*/r0*2.  The  choice  of  p*(n*)  which  matches  the 
ambient  atmosphere  at  r*"2 <  <  1  is  (upon  restoration  of  the  subscript  unity) 


pl*  =  P lmb  [z*  +  V2/(2g*  r*2)]  . 


(3.2.8) 


In  this  approximation  the  ambient  formally  is  taken  to  hold  at  an  infinite  radial 
distance,  but  the  discrepancy  from  a  match  to  an  ambient  at  a  value  of  (ro*/r-j*) 
of  only  ten  results  in  an  error  of  about  merely  one  percent. 

The  more  general  case  of  a  vortex,  the  angular  momentum  of  which  decreases 
linearly  from  a  value  of  rQ*  at  sea  level  to  the  value  of 
[(M-l )/M]rQ*,  M  ?  1 ,  at  the  tropopause  z*  =  zt*,  i.e., 

M  z  *  -  z*  M  z *  r  *2/(M  z  *) 

r*(  7*  \  a  f  *  _  n*  a  _ L  +  _2 _ P  1  "> 

1  'z  '  *0  M  z  *  n  M  zt*  -  z*  2g*  r*2 


(3.2.9) 


such  that  as  r*  -*■  “, 


z*  =  M  zt*  (n*  -  1 )/n*  . 


(3.2.10) 


Thus,  a  solution  of  the  first-order  partial  differential  equation,  p*(n*)>  any 
function  of  n*,  also  satisfies  the  boundary  condition. 


(3.2.11 ) 
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P-j *(n*)  ■*■  p*mb(z*)  for  r* 


if  one  takes  as  the  function 


P,*(n*>  =  P* 


r0*2(n  2t*  -  z*) 

2g*  r*2  M  z^* 
T0*Z(M  zt*  -  z*) 
2g*  r*2  M2  zt*2 


For  M  ■*  »,  this  case  degenerates  to  the  previous  one. 


More  generally,  if  r*(z*)  =  rQ*[(M  zt*  -  z*)/(M  zt*)]n,  n  ^  0,  then 


r  ( 1 -2n)  r  *2 

pamb  M  zt  (M  Zt  2  )  '  (M  zt*  -  z*)1"2"  2g*  r*2  (M  zt*)2" 


(3.2.12) 


p-j  *(n*)  = 


pamb  '  M  zt*  “  zt*  “  z*)  exP  “  r0*2/(2g*  r*2  M  zt*)  ■  ,  n  =  \ 


,n  t  o 


(3.2.13) 


Finally,  if  r*(z*)  =  rQ*  exp(-a*  z*),  then 


•,i*lri*>  * 


{  z*  +  (2a*)"1  Jin 


a*  r 
1  +  g fr™ 


exp(2a*  z* 
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(3.2.14) 


Thus,  deviation  of  the  isobars  (from  the  asymptotic  horizontal  planes  hold¬ 
ing  in  the  hydrostatic  ambient)  owing  to  rotation  is  known  formally  throughout 
region  I  for  any  physically  Interesting  stratification  of  the  angular  momentum, 
and  several  examples  have  been  presented  explicitly.  This  deviation  is  signifi¬ 
cant  in  ascertaining  the  pressure  differential  at  any  altitude  between  the  edge 
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of  the  eyewall  contiguous  to  the  vortex  of  region  I  and  the  edge  of  the  eyewall 
contiguous  to  the  eye  of  region  IV,  where  the  eyewall  is  envisioned  as  a  narrow 
layer.  More  explicitly,  in  region  III  the  pressure  on  the  streamsurface  con¬ 
tiguous  to  region  IV  must  match  the  known  eye  pressure  p*  (z*),  and  the  pres- 

eye 

sure  on  the  streamsurface  contiguous  to  region  I  must  match  the  pressure  in  the 
potential  vortex  at  the  eyewall  edge  p1*[R*(z*),z*]=p*amb{m*-1[m*(z*)-s*(r*)]}. 
These  statements  constitute  the  boundary  conditions  for  the  eyewall  calculation, 
to  which  attention  is  now  turned. 

The  symbol  R*(z*)  henceforth  denotes  the  radial  displacement  of  the  outer 
surface  of  the  eyewall  with  structure;  the  symbol  h*(z*)  denotes  the  transverse 
thickness  of  the  eyewall,  such  that  h*(z*)  is  measured  inward  from  R*(z*)  and 
perpendicular  to  R**(z*).  The  eyewall  is  taken  as  thin,  so  in  general  h*(z*) 

«  R*(z*);  under  this  inequality,  transversely  averaged  analyses  of  the  eyewall 
suffices. 

The  symbol  zj  denotes  the  altitude  down  to  which  hydrostatic,  adiabatic  com¬ 
pression  (of  air  whose  reference  thermodynamic  state  is  that  of  the  ambient  at 
zt*)  characterizes  Pgye(z*),  where  it  ir,  recalled  that  zt*  is  the  altitude  of  the 
tropopause  (assigned  here  to  be  the  ambient  altitude  at  which  the  pressure  and 
temperature  of  the  ambient  match  that  of  ground-level  air  which  has  risen  on  a 
moist  adiabat).  For  z*^z*>'0,  Pgye(z*)  is  given  by  hydrostatics,  and  by  moist- 
adiabatic  states  such  that  the  pressure  is  continuous  at  z*  =  zj.  If  zj  =  0,  then 
in  R*(z*)  >  r*  ^0,  the  dry-adiabatic  relation  holds  throughout  the  entire  depth 
of  the  troposphere  (for  a  moisture-free,  fully  inserted  eye).  It  is  perhaps 
worth  recalling  that  standard  thermohydrostatic  calculation  yields  p|ye(z*)  and 
zt*  from  specification  of  a  storm-prone  ambient  P£mt,(z*);  1n  the  course  of  such 
a  calculation  the  triplet  (PTOist,Tmoist’pmoist^  has  been  associated,  where 
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pamb^°^pmoist^2*^panib^2t*^’  and  where  subscript  "moist"  designates  the 
moist  adiabatic  based  on  ground-level  (or  sea-level)  ambient  conditions. 

It  is  convenient  to  introduce  the  concept  of  a  series  of  calculations, 

1  =  1,2,3...  for  known  P^^z*)  »T£mb(z*)  ,p|mb(z*) ,  and  given  r*(z*).  These 
calculations  are  based  on  different  completeness  of  "flushing"  of  moist-adiabatic 
air  from  the  eye  as  denoted  by  zj^ 1 ^ ,  i  =  1,2,3...  For  example,  one  could  study 

2*(1)=nir2t*  »  1.2,3.. .,N+1  ,  (3.2.15) 


where  i  =  1  would  constitute  a  fully  inserted  eye  flushed  of  moist  air,  and 
i  =  N  +  l  would  constitute  the  degenerate  case  of  an  eye-free  structure.  The 
theory  advanced  here  is  limited  to  cases  in  which  zj  is  a  large  fraction  of 
zt*;  i.e.,  attention  is  at  the  outset  confined  to  fully,  or  near  fully,  devel¬ 
oped  eyes. 

Under  the  cyclostrophic  and  potential -vortex  approximations,  the  peak 
swirl  speed  (of  a  two-cell  vortex  with  a  nonrotating  core)  supportable  by 
surface-level  pressure  deficit  between  eye  and  ambient,  is  estimated  by 


pa,nb<0)  - 


P*(1,(0)  * 

peye  v 


1 

T~ 


amb 


(0) 


(3.2.16) 


where  the  density  has  been  taken  as  fixed  at  its  ambient  value,  and  the  super¬ 
script  i  is  recalled  to  denote  the  completeness  of  eye  flushing.  For  i  =  1, 
the  swirl  speed  is  maximum,  and  this  value  is  used  for  nondimens ionalizati on 
below.  While  the  explicitness  of  (3.2.16)  is  convenient  as  a  tentative  approxi¬ 
mation,  a  more  accurate  formula  (that  accounts  for  the  variation  of  density 
with  radius  In  the  potential  vortex  and  that  acknowledges  that  the  peak  swirl  is 
achieved  above  the  surface  inflow  layer--rather  than  at  z*  =  G)  ultimately  is 
required  and  is  Introduced  below. 
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(3.2.17) 


The  reference  value  for  the  angular  momentum  is  taken  to  be 

r  *  =  r 
*o  rl  v 

where  r,*  is  the  (specified)  radius  at  which  the  surface  inflow  layer  sepa¬ 
rates.  If  rQ*  is  held  fixed  for  all  i,  then,  if 

rQ*0)  v*0)  =  rQ*  =>  r0*^  =  r0*/v*^  ,  (3.2. 18) 

it  follows  that 

ro*(1)  =  rj*  ;  rQ*^  >  r  ^  »  since  V*^  >  V*^  ,  i  >  1  .  (3.2.19) 

Also,  reference  values  for  the  streamwise  and  swirl  components  of  velocity 
in  the  eyewall  are  taken  as  follows: 

po*(i)  3  K  v*(i)  f  Vq*(i)  =  c  y*(i)  ,  where  0  <  C,  K  <  1.  (3.2.20,  3.2.21) 

It  may  be  noted  that 

2  2  2 

jqo*^  +  [V0>]  ^  [v^>  =»  (C2  +  K2)  <  1  .  (3.2.22,  3.2.23) 

The  equality  implies  no  kinetic-energy  losses  after  separation,  whereas  the 
inequality  implies  losses  (associated  with  vortex  breakdown).  While  formally 
a  large  range  of  values  for  C,K  may  be  explored,  in  fact  C  -  K  =  2  2  seems 
plausible. 

By  conservation  of  angular  momentum  for  the  eyewall 
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It  Is  consistent  with  nondimenslonal izatlon  introduced  in  Section  3.1  to  adopt: 
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(3.2.25) 


Hence, 


1  ){z)  v^(z)  =  z{^  ,  ij1^  =  — — 


*(i)  u  *(i) 


=  C  . 


(3.2.26) 


r  '  U* ' 
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By  conservation  of  mass  for  the  eyewall. 


[pm*o1St(2*>]  lf*U>‘2*>]  [k*1,)<2*>]  [h*(1)(z*)l  -  p‘)>(0)  q0*(i)  r/1*  !.„*« 1  >  . 

(3.2.27) 


where,  for  a  thin  layer,  it  is  reasonable  to  adopt 


ho*^  =  K  ro*  »  91ven  » 


where  rQ*  =  300  -  500  miles  for  a  typhoon,  rQ*  -  6  -  8  miles  for  a  tornado, 
and  k  =  0.0025  for  both.  If  one  defines 
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(3.2.31) 


(3.2.32) 
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In  fact,  the  value  of  <  is  not  well  set,  and  p£  g  (0)  does  not  well  charac¬ 
terize  the  density  of  the  eyewall  at  separation;  thus,  some  convenient  re¬ 
assignment  of  the  value  of  is  adopted  below  to  avoid  iteration  of 

transcendental  equations  in  establishing  initial  conditions.  However,  the 
essential  point  is  that  1S  a  known  constant.  It  may  also  be  worth  not¬ 

ing  that  the  density  of  the  moist  adiabat  P^01-st  is  tabulated  as  a  function 
of  p*  .$t,  and  thus  P^01-st  is  really  known  only  indirectly  as  a  function  of 
altitude  z*. 

Aside  from  gravitational  effects,  there  is  a  constant  head  across  all 
streamsurfaces  of  the  separating  inflow  layer  that  constitutes  the  turnaround 
and  eyewall  (region  III).  Thus,  Bernoulli's  equation  for  the  moist-adiabatic 
locus  appropriate  for  region  III  becomes  (g*  is  the  magnitude  of  the  gravita¬ 
tional  acceleration) 


[q*^(z*)~j2  +  P>(Z*)12 
L  ■  J  _  L  j  +  a*  U*  +  Z*'l  + 


p*(l)(z*) 


p*Uj(h*) 
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Jmoist  V|J2 


(3.2.33) 


where  z*’  characterizes  the  height  on  the  eye  side  of  the  separated  layer  that 
corresponds  to  height  z*  on  the  eyewall  side,  in  view  of  the  fact  that  h*  is 
measured  perpendicular  to  R**(z*)  and  that  h*  is  not  trivial  in  magnitude  for 
the  typhoon  case.  The  formula  for  z*’  is  given  below;  here  it  is  noted  that 
at  separation  z*  *  h  *^,  z*'  =  0;  also,  the  gravitational  potential  is  being 
assigned  at  midheight  of  the  layer.  As  previously  noted,  aside  from  the  gravi¬ 
tational  potential,  the  head  would  be  constant  across  virtually  all  streamsur¬ 
faces  In  the  separating  layer,  so  while  (3.2.33)  is  used  here  In  an  integral 
sense  across  the  layer,  for  the  tornado-like  case  it  would  hold  in  detail  as  well. 

The  prime  superscript  here  distinguishes  two  values  of  the  axial  coordinate,  and 
is  distinct  from  use  of  the  prime  to  denote  ordinary  derivative  in  Section  3.1. 


In  (3.2.33), 


P* 


(1) 


(z*) 


1 

~r 


P*  w  ■ 
Kamb 

m*"1 

m*(z*)  -  s*  (R*(i)(z*)) 

■  +  p*  ^ 1 ^ ( z* ' ) 
Keye  v  ‘ 

. 

. 

L—  — 

(3.2.34) 


where  p*^(z*)  could  be  written  p*^lL(z*)»  and  p*  [R*^ i ^ (z*) ,z*]  could  be 
written  p-|*[R*'^(z*),z*]. 

If  one  nondimensional izes  by  defining 
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Finally,  the  balance  of  pressure  force  and  accelerations  acting  trans¬ 
versely  to  the  separated  layer  yields 
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where 


R  (  1  )  (  7  ) 

z'  =  z  +  h^(z)  - Lt - 


{h*Pt(z)]27£ 


(3.2.39) 


Since  z  is  the  altitude  at  the  eye-contiguous  side  of  the  layer  that  is  asso¬ 
ciated  with  altitude  z  at  the  potential-vortex-contiguous  side  of  the  layer, 
the  distinction  of  z  and  z  may  be  significant  for  thick  layers  (i.e.,  for  the 
typhoon  case).  As  noted  below  (3.1.1),  the  eye-edge  pressure  exceeds  the 
ambient-edge  pressure  during  the  turnaround  (R<1),  so  the  pressure  difference 
plus  the  acceleration  owing  to  the  swirl  are  balanced  by  the  acceleration  owing 
to  the  radial-axial  flow.  The  above  formulation  is  convenient  for  the  more 
vertical  portions  of  the  turnaround,  in  which  Rz  is  zero  or  finite.  However, 
near  R=1  (i.e.,  near  separation  or  near  the  end  of  the  incursion),  |Rz|-»-°°, 
and  a  computationally  more  convenient  form  of  the  equation  is  obtained  by  inter¬ 
changing  the  roles  of  R  and  z,  such  that  the  erstwhile  dependent  variable 
becomes  the  independent  variable.  Specifically,  if  one  writes 


z  =  z 
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(3.2.40) 


then  just  after  separation,  upon  discarding  the  superscript  i  , 
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This  form  is  convenient  as  R  decreases  from  unity  until  (say)  zR  -1 ,  at 
which  point  the  form  (3.2.38)  may  be  used;  after  the  minimum  value  of  R  is 
achieved,  and  (say)  Rz  -*■  1 ,  one  may  continue  on  to  higher  z  by  use  of  the  form 


Z4  ~  P*e  (R»z)  ,  V2 

pmoist  h  (1  + 


^ZRR 


+  — - c  + - - *  (3-2-42> 

(1  +  ^  R(1  +  2^) 2  (1+z*)'2 


The  set  (3.2.26),  (3.2.30),  (3.2.36),  and  (3.2.38)  can  be  reduced  by  substi¬ 
tution  of  the  first  three  into  the  last,  to  form  one  second-order  nonlinear  ordi¬ 
nary  differential  equation  for  R(z);  admittedly,  the  factor  q(z)  is  expressed 
in  terms  of  an  intricate  functional  of  R(z),  but  the  reduction  is  still  conveni¬ 


ent  for  numerical  integration  on  a  high-speed  digital  computer.  After  speci¬ 
fication  of  two  boundary  conditions,  (say)  on  R  and  Rz  at  some  starting  value 
of  z,  the  initial-value  problem  may  be  marched  forward  in  z.  The  form  of 


(3.2.38)  convenient  at  separation  is  (3.2.41),  in  which  zR,  zRR  -*•  0  at  R  =  1 , 
z  =  h,  z'  *  0.  In  dimensional  form,  for  i  -  1  and  r*(z*)  =  rQ*,  from  (3.2.8), 


.(»)  -  ’Joist  9*  V  -  P'mb  <V  *  2^ 


(3.2.43) 


2  2  Ml2 

where  (rQ*  /r-j*  )  =  [V*v  ']  and  here  P^01-st  is  the  (known)  density  on  the  moist 
adiabat  associated  with 


ist  2  ^peye^  +  pamb  ^ho*  + 


(3.2.44) 


These  relations  give  the  sel fconsistent  means  of  assigning  V*^,  though  the 
value  varies  modestly  from  that  obtained  by  (3.2.16).  For  the  turnaround,  inte¬ 
gration  is  initiated  by  taking 


R  =  1 :  z  *  h,  zR  =  e(<0)  ,  0  <  |e|  <  <  1 , 


(3.2.45) 


or  by  using  a  locally  valid  series  approximation. 


Thus,  conservation  of  angular  momentum,  conservation  of  mass,  Bernoulli's 
equation,  and  a  balance  of  forces,  all  enforced  in  an  integral  sense,  consti¬ 
tute  four  equations  for  the  four  dependent  variables  R^(z),  v^(z),  q^(z), 
and  h^(z).  Five  dimensionless  parameters  must  be  specified:  Zj,  *  £3. 

£4,  and  further,  the  thermodynamic  and  angular-momentum  stratification 

of  the  storm-prone  (convectively  unstable)  ambient,  e.g.,  relative  humidity 
and  temperature  circulation  as  a  function  of  pressure  through  the  depth  of  the 
troposphere,  must  be  specified.  Preliminary  standard  thermohydrostatic  calcu¬ 
lations  give  zt,  Peye^(*).  Pmoist  ^pmoist^’  and  the  peak  swirl  speed  achieved 
in  the  potential  vortex  for  a  fully  inserted,  nonrotating  eye.  Specification  of 
the  five  above-cited  dimensionless  parameters  actually  entails  specification  of 
the  thickness  of  the  inflow  layer  at  separation,  h  the  radial  distance  from 

the  axis  of  symmetry  at  which  the  surface  inflow  layer  separates,  r^*;  the  height 
to  which  moist-adiabatic  air  remains  near  the  axis,  z$^;  and  the  fractions  C 
and  K  characterizing  distribution,  of  kinetic  energy  of  the  inflow  layer  at  sepa¬ 
ration,  into  swirl  speed  and  radial-axial-flow  speed. 

If,  in  the  above  one-layer  formulation,  one  has  the  same  thermodynamic  strati 
fication  and  angular-momentum  stratification  of  the  ambient,  the  same  degree  of 
flushing  of  the  eye,  the  same  ratio  of  inflow-layer  vertical  thickness  to  inflow- 
layer-separation  radius  for  a  fully  inserted  eye  (h*^VrQ*^)  and  the  same  frac¬ 
tional  distribution  (of  kinetic  energy  at  separation)  into  swirl  speed  and  radial- 
axial-flow  speed,  then  the  only  distinction  in  the  dimensionless  formulation 
between  two  cases  (one  of  which  conceivably  could  be  a  hurricane  case  and  the 
other  a  tornado  case,  within  the  above  constraints)  is  the  value  of  the  dimension¬ 
less  parameter  l ^  defined  in  (3.2.37).  The  radial  distance  to  peak  swirl  in  the 
potential  vortex,  r  =  r^*,  would  be  quite  different  in  the  two  cases;  Z3 


would  be  two  orders  of  magnitude  smaller  in  the  tornado  than  in  the  hurricane. 

This  statement  simply  presents  formally  the  earlier  observation  that  the  gravita¬ 
tional  potential  enters  significantly  only  in  the  tropical-cyclone  case.  The 
upshot  is  that  the  constant-density,  gravitational-potential-free  conditions  for  a 
periodic  solution  in  z (E^  =  0,  z  =  0  in  the  argument  of  p  ,  z  =  0  in  the  argu¬ 
ment  of  p  )  satisfactorily  hold  during  the  entire  incursion  R  <  1  and  over  some 
altitude  during  the  subsequent  excursion  R  >  1,  but  such  conditions  do  not  hold 
for  a  typhoon.  Also,  even  the  simplistic  analysis  in  Section  3.1  suggests  that 
the  separated  layer  becomes  unstable  and  undergoes  a  breakdown  during  the  excur¬ 
sion  R  >  1,  i.e. ,  at  an  altitude  of  order  r^*  or  on  the  order  of  170  m  off  the 
ground  for  a  tornadic  mesocyclone.  However,  variable-density  and  gravitation- 
potential  effects  enter  significantly  over  the  appreciably  greater  axial  dis¬ 
tance  pertinent  for  the  turnaround  in  a  typhoon  (here  rj*  *  10-16  km),  and 
entirely  different  results  are  obtained.  The  results  indicate  that  the  eyewall 
thickness  in  a  typhoon  increases  so  much  at  midtropospheric  heights  that  the 
postulated-thin-layer  basis  of  the  analysis  (R»h)  is  seriously  violated;  the 
model  then  is  inadequate. 

3.3  Numerical  Examples 

In  support  of  the  statements  at  the  end  of  Section  3.2,  calculationsare 
reported  that  are  based  on  the  mean  ambient  for  the  West  Indies  for  September 
(i.e.,  for  the  peak  of  the  hurricane  season)  (Jordan  1957).  Figure  6  presents  the 
temperature  as  a  function  of  pressure  for  ambient  air,  for  sea-level  ambient  air 
lifted  on  a  moist  adiabat  (taken  to  be  an  idealized  characterization  of  the  thermo¬ 
dynamic  states  of  the  eyewall),  and  for  tropopause-level  air  dry-adiabatical ly 
compressed  down  to  sea-level  altitude  (taken  to  be  an  idealized  characterization 


32 


of  the  thermodynamic  states  of  the  eye);  constant-altitude  lines  are  assigned  on 
the  basis  of  hydrostatics.  From  these  thermohydrostatic  calculations. 


zt*«13.8  km,  P*mb(0)  =  1014  mb,  p|ye(0)  =  91 7  mb,  P*ye(0)  =  9.35 x  10‘4  g/cm  ; 
from  (3.2.16) 


V*^  =  131.6  m/s. 


(3.3.1a) 

(3.3.1b) 


Plausible  assignment  of  lateral  scales  suggests 

rx*  =  1.52  x  104  m  ,  <  =  2.5  x  10"3,  rQ*  =  6.09  x  105  m.  (3.3.2) 
Thus,  from  (3.2.18)  and  (3.2.28), 

rQ*=  rx*  V*(1)  =  2.01  x  106  m2/s,  hQ*(l)  =  k  rQ*=  1.72x  103  m  .  (3.3.3) 

From  (3.2.20)  and  (3.2.21),  if  C  =  K  =  2"^,  i.e.,  for  no  kinetic-energy  loss, 
qo*(D  s  k  v*(1)  =  92.9  m/s  ,  vQ*(1)  «  C  V*(1)  =  92.9  m/s  (3.3.4a) 

Alternatively  (and  less  plausibly),  for  about  25%  loss  of  kinetic  energy  prior 
to  any  vortex  breakdown,  if  K  =  0.7  and  C  =  0.5, 

qo*U)  =  91.4  ro/s  ,  Vq*(1)  =  65.3  m/s  .  (3.3.4b) 

Accordingly,  under  (3.3.4a),  since  r^*  =  r  by  definition, 

Z^1*  *  0.707,  Z2^  =  8.83  x  10-2 ,  Z3  *  17.5,  Z4  =  11.4,  Z^1*  =  1,(3. 3. 5a) 

where  these  dimensionless  groups  are  defined  by  (3.2.26),  (3.2.32),  and  (3.2.37) 
Under  (3.3.4b), 

=  0.5,  Z =  7.49  x  10" 2 ,  Z,  =  17.5,  Z,  =  11.4,  I *  0.74. 


Under  (3.3.4b)  and  (3.3.5b),  for  a  fully  inserted  eye,  the  results  of  numeri¬ 
cal  integration  of  the  initial-value  problem  for  the  turnaround  (described  in  the 
penultimate  paragraph  of  Section  3.2)  is  given  in  figures  7  and  8.  Actually,  as 
discussed  in  the  last  paragraph  of  Section  3.2,  some  restrictive  parametric 
assignments  are  made  to  permit  simulation  of  the  tornado:  2^=0,  z  =  0  in  the 
argument  of  p  (so  p  =  1).  z  =  0  in  the  argument  of  p  as  given  by  (3.2.8), 
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and  z  =  0  in  (3.2.45).  (The  hurricane  ambient  is  seen  below  not  to  be  very  dif¬ 
ferent  from  some  tornado  ambients.)  Under  the  simplifications  just  enumerated, 

the  compressible  formulation  of  the  turnaround  (developed  in  Section  3.2)  does 
degenerate  essentially  to  the  incompressible  treatment  (developed  in  Section 
3.1).  Explicitly,  there  is  an  incursion  (from  R  =  1)  of  the  separated  layer 
toward  the  axis  of  rotation  such  that  there  is  a  moderate  overshoot  of  the 
swirl,  then  an  excursion  such  that  an  instability  leading  to  vortex  breakdown 
occurs  (figure  7),  all  over  an  altitude  of  z  =  0(1);  furthermore,  the  eye-edge 
pressure  exceeds  the  pressure  at  the  ambient-edge  pressure  during  incursion  (the 
opposite  being  true  during  excursion),  with  the  thickness  h  and  radial-axial- 
flow  speed  q  being  rather  constant  with  height  (figure  8). 

Under  (3.3.4a)  and  (3.3.5a),  again  for  a  fully  inserted  eye,  but  without, 
any  of  the  restrictive  assignments  listed  in  the  last  paragraph  to  permit  simu¬ 
lation  of  the  tornado,  one  finds  results  holding  for  the  typhoon  (figures  9-11). 
Particularly  noteworthy  is  the  very  rapid  increase  in  thickness  for  z  =  0.35, 
such  that  (h/R)  =  0.5  and  is  increasing  rapidly  with  altitude;  the  bases  for  the 
integral  treatment  have  deteriorates 


4.  The  Eyewall 
4.1  General  Comments 

The  Integral  theory  for  the  evolution  of  the  separated  layer  (region  III), 
after  the  postulated  breakdown  of  the  turnaround  discussed  in  Section  3.1,  is 
now  given  for  the  tornadic  mesocyclone.  The  four  basic  equations  remain  those 
presented  earlier:  (3.2.26),  (3.2.30),  (3.2.36),  and  (3.2.38).  However,  while 
the  independent  variable  remains  the  distance  above  ground  level  z,  it  has 
already  been  anticipated  (in  the  last  paragraph  of  Section  2.2)  that  the  roles 
of  certain  specified  and  inferred  variables  would  be  interchanged.  Specifically, 
only  certain  axial  profiles  of  peripheral  angular  momentum  are  analytically 
tractable  for  closed-form  presentation  in  the  characteristic  argument  of 
p(i)(R(i),z) ^  and  study  of  these  particular  profiles  often  leads  to  results  incon- 
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sistent  with  two-cell  structure.  For  example,  R^(z)  may  become  large  and  nega¬ 
tive  at  midtropospheric  heights,  upon  integration  forward  in  z  after  adoption  of 
suitable  starting  values  (see  below);  such  behavior  (deflection  of  the  eyewall 
toward  the  axis)  is  Inconsistent  with  eye  insertion,  since  the  eye  extends  down¬ 
ward  from  the  neutral -^stabll  ity  altitude.  It  turns  out  to  be  much  more  compu¬ 
tationally  efficient  to  specify  R^(z),  the  displacement  of  the  separated  layer 
from  the  axis  of  rotation,  and  to  use  the  four  equations  to  obtain  pfl^(R^Hz), 
q^(z),  v^(z),  and  h^(z)  for  a  given  ambient  thermodynamic  stratification 
[such  that  pgyg(z)  is  given].  In  a  separate,  sequential  calculation  one  may 
infer  the  associated  r^(z)  by  inversion  of  a  triangular  matrix  (see  below); 
in  this  manner,  one  can  obtain  a  (piecewise-constant)  profile  for  r^(z)  for 
those  cases  for  which  two-cell  structure  is  consistent  with  the  model. 


From  inspection  of  (3.2.26),  (3.2.30),  (3.2.36),  and  (3.2.38),  and  from 
plausible  physical  behavior  for  two-cell  structure,  one  anticipates  that  R^(z) 
increases  monotonically  in  z,  where  R^(z)  >  1;  that  p  ^  t(z),  v^(z),  and 
p^(z)  decrease  monotonically  in  z,  where  p^(z)  =  Cpe^e(z ' )  +  p^^(R^\z)]/2; 
and  that  q^(z)  and  h^(z)  may  either  increase  or  decrease  with  z,  and  there¬ 
fore  may  not  be  monotonic  in  z.  Since  h^(z)  is  small  for  the  tornadic  case, 
distinction  of  z  and  z  seems  less  critical  for  this  case. 

4.2  Computational  Details;  Starting  Conditions 

The  partially  algebraic,  partially  differential  set  conveniently  is  reduced 
to  a  purely  algebraic  set  upon  specification  of  R^(z).  Nevertheless,  at  each 
height  z,  a  double  iteration  on  h^(z)  and  p^(z)  is  required  because  of  the 
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nonlinear,  coupled  nature  of  the  equations.  In  particular,  upon  specification 
of  R^(z),  from  which  v^(z)  follows  immediately  by  (3.2.26),  adoption  of 
trial  values  for  h^(z)  and  p^(z)  permits  obtaining  a  compatible  value  for 
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q^(z)  from  (3.2.36).  In  this  connection,  the  inputs  z  and  p ^(z)  are 
readily  available  from  (3.2.39)  and  from  tabulated  results  from  standard  thermo¬ 
hydrostatic  calculations.  Updated  values  for  pf^(z)  and  for  h^(z)  are  then 
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sequentially  available  from  (3.2.38)  and  (3.2.30),  respectively.  One  than  may 
update  q^(z).  Effective  invariance  between  successive  iterates  for  all 
variables  is  taken  to  indicate  convergence;  the  value  of  z  is  incremented,  and 
the  procedure  repeated.  It  may  be  remarked  that  adopting  straight-line  seg¬ 
ments  for  R^(z),  with  neglect  of  the  right-hand  side  of  (3.2.38)  is  a  plausible 
procedure  since  one  anticipates  very  modest  curvature  over  most  of  the  eyewall; 
substitution  of  a  smoothed  profile  for  R^(z)  should  then  be  carried  out  to 
confirm  the  relative  unimportance  of  the  omitted  term. 
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The  boundary  conditions  are  taken  to  be  as  follows: 


z  =  0:  R(i)=  1,  Rz(l)=  r£).  0  .  (4.2.1) 

These  conditions  imply:  (1)  the  altitude  over  which  the  turnaround  occurs  for 
the  tornadic  case  is  so  much  less  than  the  height  of  the  neutral -stability  lid 
that  the  eyewall  might  just  as  well  be  regarded  as  starting  from  ground  level 
z  =  0;  (2)  the  vortex  breakdown  occurs  so  soon  after  the  completion  of  the  turn¬ 
around  that  setting  R^=  1,  i.e.,  starting  the  eyewall  at  the  radial  distance 
of  inflow-layer  separation,  is  a  plausible  procedure;  and  (3)  the  separated  layer 
emerges  from  the  breakdown  in  nearly  vertical  ascent. 

A  significant  constraint  on  physically  acceptable  profiles  for  R^(z),  not 
previously  discussed,  is  the  compatibility  with  the  stability  criterion  in  swirl¬ 
ing  flows  that  (at  each  altitude)  the  angular  momentum  not  decrease  with  increas¬ 
ing  radial  distance  from  the  axis  of  rotation.  In  particular,  since  in  the  inte¬ 
gral  approach,  the  transversely  averaged  value  v^(z)  characterizes  the  eyewall 
rotation  at  each  height,  at  the  interface  r=R  between  the  eyewall  (region  III) 
and  the  potential  vortex  (region  I),  it  is  required  that  v^(z)  <  T^(z)/R^(z). 
That  is,  the  displacement  of  the  eyewall  annulus  from  the  axis  of  symmetry  must 
be  such  that  the  swirl  of  the  fluid  in  the  eyewall  is  less  than  the  swirl  of  the 
fluid  contiguous  to  the  eyewall  in  the  potential  vortex.  Only  high  in  the  tropo¬ 
sphere,  at  altitudes  at  which  outflow  occurs,  may  this  inequality  be  violated. 
Although  discussion  is  deferred  to  Section  4.3  concerning  how  one  obtains  the 
profile  of  ambient  circulation  with  altitude,  r^^(z),  that  is  compatible  with 
the  eyewall  standoff  R^(z)  adopted,  nevertheless  this  important  constraint  on 
R^(z)  should  be  noted. 


It  is  of  physical  interest  to  obtain  solutions  for  the  following  parametric 
variations:  for  several  measured  profiles  of  temperature  and  relative  humidity 
with  pressure  for  tornado-prone  ambients;  for  several  degrees  of  eye  insertion  in 
each  such  ambient,  i  =  1,2...;  and  (perhaps)  for  several  characterizations  of 
the  kinetic-energy  loss  through  the  breakdown,  K,C,  though  it  may  be  that  K  =  C  *  2^ 
is  the  case  of  dominant  physical  interest. 

It  may  be  apropos  to  note  that  the  conditions  for  the  radial  displacement 
of  edge  of  the  eyewall  contiguous  to  the  eye,  given  by 

R(i)(z)  -  h(i)(z)/jl  +  R*1}(z)  2|2  ,  (4.2.2) 

to  be  constant  in  z,  is  of  physical  (and  historic)  interest.  The  present  model 
loses  self consistency  as  this  condition  is  approached,  because  the  augmented 
pressure  deficit  from  ambient  owing  to  compressional  heating  of  the  air  contigu¬ 
ous  to  the  axis  would  become  unavailable  to  sustain  swirling  at  speeds  in  excess 
of  those  compatible  with  moist-adiabatic  ascent. 

Finally,  for  the  tornado-eyewal 1  calculation,  assignment  of  the  peak  swirl 
V*^1)  by  use  of  (3.2.43)  incurs  a  small  error,  which  in  turn  leads  to  a  minor 
inconsistency  in  the  angular-momentum  field.  As  noted  in  the  second  paragraph 
of  this  section,  the  turnaround  is  conveniently  ignored,  and  the  eyewall  is 
taken  to  rise  vertically,  at  a  prescribed  distance  from  the  axis  of  symmetry, 
at  the  starting  altitude  z  *  0.  Thus,  the  assignment  of  V*^  at  the  inflow- 
layer-separation,  as  carried  out  in  (3.2.43),  is  inappropriate.  Indeed,  compari¬ 
son  of  (3.2.38),  which  is  taken  to  hold  at  z'  =  z  =  0,  with  (3.2.43)  is  informa¬ 
tive.  Since  Rz  »  RZ2  =  0  and  the  term  involving  Z3  is  negligible,  as  is  the  gravi¬ 
tational  term  in  (3.2.43),  one  sees  that  (3.2.38)  implies  a  value  of  the  pressure 


at  the  potential-vortex  edge  of  the  thin  layer,  p  ,  that  is  greater  after  turn- 
around  than  is  the  analogous  pressure  at  inflow-layer  separation;  the  last  term 
on  the  left-hand  side  of  (3.2.38)  is  the  basis  of  the  discrepancy--swirl  is  a 
source  of  pressure  difference  across  the  post-turnaround  thin  layer,  but  not 
across  the  thin  layer  at  separation.  The  upshot  is  that  V*^  is  to  be  found  by 
iteration  by  the  requirement  that  the  dimensionless  angular  momentum  achieve 
its  nominal  value  of  unity  at  the  (pseudo-)  ground-level  position  z  =  0.  Of 
course,  (3.2.43)  provides  a  good  first  estimate  for  iteration.  At  this  point, 
discussion  of  how  the  angular-momentum  profile  is  obtained  is  appropriate. 


4.3  Angular-Momentum  Distribution  with  Altitude  Consistent  with  a  Prescribed 

Eyewall  for  a  Thermodynamically  Specified  Ambient 

Via  the  characteristic  curve  delineating  the  locus  of  an  isobar  extending 
from  the  outer-eyewall  edge  to  the  ambient  in  the  hydrostatic,  cyclo- 
strophic,  potential  vortex  (region  I),  one  can  complete  the  inverse  solutions; 
i.e.,  one  can  identify  the  peripheral  circulation  compatible  with  an  eyewall 
under  prescribed  axial  profiles  for  radial  displacement  from  the  axis,  thickness, 
swirl,  radial-axial  flow,  potential -vortex-edge  pressure,  etc. 

Since  the  identification  is  to  be  numerical,  the  sought  function  x(z), 
where  (recall  rQ*  *  rQ*^  V*^) 


x(z)  =  **  [x(z)]' 


(4.3.1) 


is  found  in  piecewise-constant  approximation.  Explicitly,  if  H(z)  denotes  the 

Heaviside  unit-step  function, 

M 

x(z)  -  l  xk  [H(z  -  zk)  -  H(zk+1-  z)]  ,  xk  const., 
k— 1 


(4.3.2) 


where  for  convenience 


zk  *  zM  -  (M-k)(Az),  k  =  M,(M-l),(M-2) . 1  , 


(4.3.3) 


so  zk+1  >  zk,  Zj  >  0  and  M  is  the  number  of  discrete  values  used  to  sample  the 
eyewall  solution: 


»  (p^g)  >  hk  ,  qk  ,  vk  ,  ...  vs  zk  ,  k  —  1,2,3,...,M  .  (4,3.4) 


[In  this  section  the  superscript  i,  introduced  previously  to  distinguish  dif¬ 
ferent  degrees  of  eye  insertion  for  a  given  thermodynamic  stratification  of 
the  ambient,  is  taken  as  implicit,  for  simplicity  of  notation.]  It  turns  out 
that  choosing  an  interval  (Az)  and  inferring  M  is  preferred.  If  ZM  is  defined 
by 


^amb^M^  -  Pae^zM^  ’ 

where  zM(<zt)  is  the  greatest  height  to  which  the  eyewall  solution  is  carried 
(or  is  given  creditability),  and  if  N  .  is  the  minimum  positive  integer  satis¬ 
fying  the  inequality 


(Nmin)(Az)  >  (ZM  *  ZM> 


(4.3.5) 


then  it  is  taken  that 

42  2  <\1nK42)  . 


(3.4.6) 


Then,  from  (4.3.3),  zk  may  be  assigned,  where  z-j  is  the  smallest  positive  value 

assigned.  The  values  cited  in  (4.3.4)  may  be  found  by  interpolation  in  the 

tabulated  values  for  the  eyewall  solution;  in  particular,  the  values  (p  )  are 

ae  k 

pertinent  since  one  needs  to  tabulate 


lWZk>  ■  Pae<zk>  *  (P«>,  •  K  = 


Another  quantity  required  below  is 


Int 


zk  -  zk 


Az 


=  Nk  ,  k  =  1,2,3,. .. ,M  , 


(4.3.7) 


(4.3.8) 


where  the  symbol  "Int"  indicates  the  largest  positive  integer  less  than  or 
equal  to  the  argument.  In  general,  N"k  is  nonincreasing  in  value  with 
increasing  k. 


In  effect,  Gaussian  elimination  (back  substitution)  is  being  used  to  solve  a 
set  of  simultaneous  linear  algebraic  equations  with  a  triangular  coefficient 
matrix.  It  is  recalled  that  rk  =  s,  where  zk+1  -  =  (Az)  and  zM,  the  peak 

value  of  z,  is  the  highest  value  utilized  from  the  eyewall  solution. 

It  seems  worth  recalling  that  in  (4.3.11)  the  position  z  =  0  corresponds 
to  the  height  at  which  the  turnaround  is  completed.  It  is  remarked  that  no 
solution  is  being  furnished  for  the  angular-momentum  profile  between  ground 
level  in  the  ambient  and  that  height  in  the  ambient  at  which  (essentially)  the 
pressure  has  decreased  to  the  pressure  at  the  base  of  the  eye. 


a 


m 


4.4  Numerical  Examples 
4.4.1  Strong  Tornado 

One  of  the  most  convectively  unstable  ambients  ever  measured  in  the 
vicinity  of  a  tornado  was  taken  near  St.  Cloud,  MN,  on  July  22,  1967  at  1715Z 
(figure  12).  At  ground  level,  Ta*b(0)  =  307.2  K,  Pa*5(0)  =  1.12x10* 3  g/cm3, 
Pa*b(0)  =  1002.4  mb.  The  pressure  at  the  base  of  a  moist  adiabat  (idealized 
core  of  a  one-cell  vortex)  is  calculated  to  be  930.4  mb,  such  that  the  peak 
swirl  for  a  Rankine  vortex  sustained  by  this  pressure  deficit  would  be  80.0  m/s. 
The  pressure  at  the  base  of  a  fully  inserted  eye,  idealized  as  dry-adiabatical ly 
compressed,  is  calculated  to  be  P  *  ^(0)  s  839.7  mb,  such  that  the  peak  swirl 

CJT  C 

for  a  potential  vortex  V*^  =  170.3  m/s.  The  tropopause  height  z£  =  13.2  km. 
The  thermohydrostatic  data  for  this  ambient  are  presented  in  figure  13.  It 
should  be  emphasized  that  the  idealized  processes  adopted,  together  with  this 
singularly  unstable  ambient,  lead  to  exceptionally  large  peak  speeds  and  maxi¬ 
mum  pressure  deficits. 


The  following  parametric  assignments  are  adopted: 

r*  =  152  m,  h  *  =  25.4  m,  K  =  C  =  2~h. 


(4.4.1) 


Iteration  reveals  that  the  selfconsi stent  value  for  V*'  discussed  in  the  last 


paragraph  of  Section  4.2  is  given  by 

u*(l)  . 


164.0  m/s  , 


(4.4.2) 


or  about  3.5%  below  the  value  obtained  by  applying  cyclostrophy  to  the  thermo- 
hydrostatically  computed  pressure  deficit  at  the  base  of  the  eye.  This  slightly 
reduced  value  of  V*^  yields  r(0)  =  1.  If  the  superscripts  i  are  discarded, 
it  follows  that  the  derived  nondimensional  groups  take  on  the  following  values: 


V  V  V  V  V  *-  ' 


as 

< 

xx 


zl  =  0.707,  Z2  =  0.141,  Z3  =  0.111,  Z4  =  7.55,  Zg  =  1.0  .  (4.4.3) 

For  a  fully  inserted  eye,  firgure  14  presents  the  eyewall  geometry  as  a 
function  of  altitude  z  above  ground  level  (more  meticulously,  above  the  height 
for  completion  of  the  turnaround).  The  outer  eyewall  locus  R(z)  is  found,  by 
trial  and  error  (see  below),  in  the  form  of  straight-line  segments,  and  the 
inner  eyewall  is  then  deduced  by  (4.2.5).  It  has  been  confirmed  that  smoothing 
the  segmented  profile  gives  very  nearly  the  same  result;  i.e.,  the  contribution 
from  the  term  on  the  right-hand  side  of  (3.2.38)  is  negligible  relative  to  the 
largest  terms  on  the  left-hand  side.  Figure  15  presents,  as  a  function  of 
height  z,  the  eyewall  thickness  h  and  radial-axial-flow  speed  q,  for  the  adopted 
function  R(z);  these  quantities  vary  modestly  with  height,  and  resotration  of 
eye-eyewall  entrainment  may  be  anticipated  to  counteract  the  depicted  thinning 
of  the  eyewall  in  the  upper  troposphere.  Figure  16  presents  the  merging  of  the 
pressure  in  the  eye  and  the  pressure  at  the  ambient  edge  of  the  eyewall,  as  the 
potential  isobaric  impervious  "lid"  at  the  tropopause  is  approached.  Figure  17 
presents  the  distribution  of  the  angular  momentum  r  as  a  function  of  altitude  z; 
the  key  observation  here  is  the  constancy  of  the  antular  momentum  at  the  (appre¬ 
ciable)  lower- tropospheric  value  through  much  of  the  troposphere.  Figure  17  pre¬ 
sents  the  transversely  averaged  eyewall  swirl  speed  v,  and  the  swirl  speed  in 
the  potential  vortex  at  the  interface  between  the  potential  vortex  and  the  eye- 
wall  r/R,  as  functions  of  height  z.  For  stability  it  is  necessary  that  (r/R)>v, 
and  it  is  seen  that  this  inequality  is  violated  only  very  high  in  the  tropo¬ 


sphere.  It  is  the  satisfaction  of  this  inequality  that  necessitates  lengthy 
trial-and-error  numerical  experimentation  with  profiles  for  R(z). 


For  an  eye  that  is  only  90%  inserted,  i.e.,  for  the  case  in  which  (for  the 
same  ambient)  the  bottom  10%  in  height  of  the  eye  consists  of  moist-adiabatic 


air  (such  that  the  pressure,  but  not  density  or  temperature  or  mixing  ratio, 
is  continuous  at  the  dry-air/moist-air  interface),  analogous  plots  of  results 
are  presented  in  figures  18-21.  The  alterations,  mainly  at  low  altitude,  are 
selfevident  and  modest;  consequences  of  the  density-discontinuity  interface  in 
the  eye  appear  in  some  of  the  profiles.  For  this  partial  insertion,  the  peak 
swirl  is  reduced  to  155.8  m/s;  for  invariance  of  the  angular  momentum  rQ*  = 
rj*  V*^,  the  radius  of  peak  swirl  increases  to  160.4  m.  Whereas  Zj,  E^,  and 
are  unaltered  from  the  values  given  in  the  last  paragraph  for  the  fully  in¬ 
serted  case,  ^2  ~  0.144,  =  0.902.  The  principal  point  is  that  the  theory 

does  seem  robust  enough  at  least  to  begin  inquiry  about  less-than-complete  eye 
insertion. 

4.4.2  Moderate  Tornado 

A  moderately  convectively  unstable  ambient  measured  in  the  vicinity  of  a 
tornado  was  taken  near  Jackson,  MS,  on  April  17,  1978  at  2300Z  (figure  23).  At 
ground  level,  T*mb(0)  =  301.2  K,  p*mb(0)  =  1.158xl0-3  g/cm3,  P*mb(0)  =  1009.9 
mb.  The  pressure  at  the  base  of  a  moist  adiabat  (idealized  core  of  a  one-cell 
vortex)  is  calculated  to  be  981.2  mb,  such  that  the  peak  swirl  for  a  Rankine  vor¬ 
tex  sustained  by  this  pressure  deficit  would  be  49.6  m/s.  The  pressure  at  the 
base  of  a  fully  inserted  eye,  idealized  as  dry-adiabatically  compressed,  is  cal¬ 
culated  to  be  p*  (0)  =  923.9  mb,  such  that  the  peak  swirl  for  a  potential  vortex 
eye 

V*U)  =  122.0  m/s.  The  tropopause  height  z*  -  12.0  km.  The  thermohydrostatic 
data  for  this  ambient  are  presented  in  figure  24.  The  less  severe  nature  of  the 
convective  instability  for  this  ambient,  re' .five  to  that  introduced  in  Section 
4.4.1,  is  evident. 

The  parametric  assignments  of  (4.4.1)  are  adopted  again.  Iteration  reveals 
that  the  selfconsistent  value  for  V*^  discussed  in  the  last  paragraph  of 


Section  4.2  is  given  by 


V*^  =  164.0  m/s  ,  (4.4.4) 

or  about  3.9%  below  the  value  obtained  by  applying  cyclostrophy  to  the  thermo- 
hydrostatically  computed  pressure  deficit  at  the  base  of  the  eye.  If  the  super¬ 
scripts  i  are  discarded,  it  follows  that  the  derived  nondimens ional  groups  take 
on  the  following  values: 

Zj  =  0.707,  Z2  =  0.133,  Z3  =  0.226,  Z4  =  14.4,  Z5  =  1.0  .  (4.4.5) 

Results  for  a  fully  inserted  eye  for  the  moderately  convectively  unstable 
ambient  of  figures  23  and  24  are  presented  in  figures  25-29;  these  may  be  juxta¬ 
posed  with  their  counterparts  (figures  14-18)  for  a  fully  inserted  eye  for  the 
highly  convectively  unstable  ambient  presented  in  figures  13  and  14.  The  eye- 
wall  is  appreciably  more  vertical  (and  somewhat  thinner)  for  the  more  unstable 
ambient.  However,  the  key  result  that  for  two-cell  structure  the  ambient  cir¬ 
culation  remains  rather  invariant  with  height,  and  retains  its  ground-level 
value  until  high  in  the  troposphere,  holds  for  the  moderately  convectively 
unstable  ambient  as  well  as  for  the  highly  convectively  unstable  ambient. 

Although  for  brevity  no  figures  are  presented  for  the  case  of  90%  eye 
insertion  for  the  moderately  convectively  unstable  ambient,  in  fact  these  calcu¬ 
lations  have  been  carried  out  and  give  the  same  trends  as  those  reported  in  Sec¬ 
tion  4.4.1  for  the  si ightly-less-than-completely-inserted  eye. 


5.  Remarks  on  Possible  Further  Directions 

Future  work  on  the  model  under  discussion  might  go  on  to  consider  the 
upper-level  outflow  of  "processed  throughput".  However,  probably  further 
refined  treatment  of  the  two-cell  structure  of  the  vortex  core  is  warranted. 

It  has  already  been  noted  that  addition  of  a  model  of  entrainment  to  the 
inviscid  analysis  undertaken  here  would  be  of  high  priority.  While  entrain¬ 
ment  models  for  rapidly  rotating  flows  are  not  very  advanced,  it  seems  likely 
that  the  swirling  and  rising  air  of  the  lower  eyewall  annulus  would  entrain 
fluid  from  the  virtually  nonrotating  eye.  Indeed,  such  entrainment  may  well 
be  necessary  for  selfconsistency  of  the  model  under  closer  scrutiny,  for  air 
from  the  virtually  nonrotating  thin  diffusive  sublayer  of  the  low-level  inflow 
region  (II)  is  continuously  being  conveyed  toward  the  axis  for  an  incompletely 
inserted  eye.  This  air  would  accumulate  in  time  and  obviate  the  postulated 
quasisteady  nature  of  the  mature-stage  two-cell  structure;  thus,  this  steadily 
supplied  low-momentum  sublayer  air  must  be  conveyed  away  by  entrainment  into 
the  eyewall . 

Even  before  entrainment  is  introduced,  it  would  seem  worthwhile  to  return 
to  the  direct  problem  of  specifying  axial  distributions  of  angular  momentum 
compatible  with  two-cell  structure.  Experience  with  inverse  solutions  may 
render  such  an  approach  feasible.  At  least,  a  perturbational  analysis  about  a 
distribution  r(z)  compatible  with  two-cell  structure  might  suggest  the  sensi¬ 
tivity  of  results  to  profile  properties. 

Finally,  the  current  modeling  of  the  eyewall  by  a  moist-adiabatic  locus  is 
overidealized.  It  seems  unlikely  that  the  ambient  stratification  of  total  static 


enthalpy  holds  uniformly  in  the  potential  vortex  (region  I)  and  the  moist 
adiabat  holds  uniformly  in  the  eyewall  (III).  Provision  for  some  smoothing 
of  this  transition  may  warrant  consideration. 


APPENDIX  A.  A  DERIVATION  OF  THE  MOIST  ADIABAT 


Consider  a  volume  V  with  surface  S  and  surface  element  dS;  the  outward 
normal  vector  to  dS  is  denoted  n^  and  fluid  velocity  vector  is  denoted  u . . 
Super  asterisks,  used  in  the  main  text  to  denote  dimensional  quantities, 
are  dispensed  with  for  this  appendix. 

For  conservation  of  mass,  consider  air  first.  The  flux  of  dry  air 
through  dS  is  in  Cartesian  notation,  p  u^  n..  dS,  where  p  is  air  density. 
Hence,  in  the  steady  state, 

p  ui  ni  dS  =  0  =>(p  Uj)j  =  0,  (A.l ) 

where  comma  denotes  partial  differentiation.  Water  vapor  has  the  same 
speed,  but  density  a.  When  a  <  E(T),  where  Z(T),  the  saturation  vapor 
density,  is  a  known  function  of  the  temperature  T, 

(a  u.j )  i  =  0. 

If  the  mixing  ratio  for  water  vapor  Y  =  o/p,  then 

p  u.  Y  .  =  0.  (A. 2) 

*  f  * 

< 

When  T  is  low  enough. 


o  =  E(T) .  (A. 3) 

For  conservation  of  momentum,  the  momentum  flux  associated  with  u^ 

is  (p  +  o)u.  n.  u.  dS,  and  that  associated  with  the  partial  pressure  of 
J  J  ■ 

air  p  and  with  the  partial  pressure  of  vapor  a  Is  (a  +  p)n^  dS.  The  pro¬ 
duction  of  momentum  per  unit  volume  associated  with  the  gravitational 
body  force  (g  Is  the  magnitude  of  the  acceleration  of  gravity)  is 

-g  z  .(p  +  o),  where  the  spatial  coordinate  z  lies  parallel  to  the  body 
#  ■ 

force.  The  loss  of  momentum  per  unit  volume  associated  with  the  dis¬ 
appearing  fluid  is  -(a  u.)  .  u^ .  Thus,  using  (A.l), 


(b  (p  +  a)u.  n.  u.  dS  +  ®  (a  +  p)n.  dS  =  -g  /  (p  +  a) z  . 
/c  J  J  '  *  •'II  >  ■ 


J  (o  u.)  .  u.  dV  (A. 4) 
•'ll  0  »3  1 


becomes 


(p  +  a)u .  u.  +  (a  +  p)  .  +  (p  +  o)g  z  .  =  0. 

J  1  :  9  •  »  • 

»  J 


(A. 5) 


Incidentally,  if  a  momentum  equation  had  been  written  for  each  constituent, 
dry  air  and  water  vapor,  interparticle  drag  of  the  form  D  =  A(ui  -  v^) 
would  have  been  introduced,  where  u^  is  air  velocity  and  v^  is  vapor  velo¬ 
city,  and  here  A  °°  so  vi  =  u^. 

The  equations  of  state  are 


p  =  P  R  T, 


(A. 6) 


a  =  a  R'  T, 


(A. 7) 


where  the  gas  constant  for  vapor  R'  is  related  to  the  gas  constant  for 
vapor  by  R'  =  R/(mv/ma),  where  my  is  the  molecular  weight  of  water  vapor 
and  is  the  molecular  weight  of  dry  air. 

The  nonmechanical  energy  carried  by  the  vapor  can  be  written  as 
| L0  +  Cy  t),  where  L0  is  the  (extrapolated)  value  of  the  heat  of  condensa¬ 
tion  at  T  =  0,  and  c^  is  the  specific  heat  at  constant  volume  of  the 
vapor  v.  With  this  description,  the  net  heat  released  by  condensation 
at  temperature  T  would  be 


L  *  L0  +  (c;  -  c)T , 


(A. 8) 


where  c  is  the  specific  heat  of  liquid  water.  The  flux  of  energy  across 
dS  via  u.|  is 


cl  O  c* 


.1  (.'  tUV, 


(A. 9) 


[»  -f  («,  T  *  ^-s-se.)  +  a  u,  (<  T  ♦  ♦  L„)J  »,  dS, 


where  cy  is  the  specific  heat  at  constant  volume  of  dry  air.  The  stream- 
wise  velocity  component  is  denoted  q  and  the  perpendicular  velocity  com¬ 
ponent  (swirl)  is  denoted  v.  Below,  Q2  h  (v2  +  q2 ) .  The  work  done  by 
pressure  forces  across  dS  is  (a  +  p)u^  n^  dS,  and  the  flux  of  potential 
energy  is  (p  +  a)g  z  n^  u..  dS.  The  rate  at  which  energy  goes  into  a  sink 
(i.e.,  into  liquid,  which,  by  conventional  approximation  1:£:reat1ng  the 
moist  adiabat,  just  vanishes)  is  the  product  of  two  factors:  -(a  u^)  ^ , 
the  rate  of  loss  of  mass  per  unit  volume,  and  jg  z  +  £(q2  +  v2 ) /2 J  +  c  t|. 
It  follows  that 

[p  uj(Cv  T  +  ¥")  +  0  uj(cv  T  +  Lo  +  +  (p  +  a)uj  3  2  +  •  +  p)ujj 


Hence, 


*  (o  ui)  T  +  |  M  2 j  . 

“j  (9  2  +  ^),j +  (a  uj!.j  [(<  -  C)T  +  Lo] 

+  0  uj(<  T  +  I-  *  9  ^  +  p  “j<c,  T>,j 


(“j  P.j  +  P  “j.j)  *  (UJ 


(A. 10) 


(A. 11) 


where  (a  uj)  j  *  p  uj  Y  ^ ,  a  u^  *  p  Y  u^,  and  the  last  line  of  (A. 11)  is 


p  uj(p),j  +  p 


(A. 12) 
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(A. 13) 


Thus,  since  (cy  +  R)  =  cp  and  (c^  +  R')  =  c^. 


Uj(l  ♦  V)  ( i:4)  *,  Uj(cp  T)  |j  +  p  Uj  V#j(Lo  -  c  T) 


+  p  u  (Y  cl  T)  .  =  0, 

J  r  t  J 


where  cv  (not  cv)  appears  because  work  is  done  on  condensing  gas  and  the 
p  *  ri  -  (c  -  OtT  But 

_  o  D  j 


enthalpy  release  is 


P  u. 


r  1  (P  +  ct)  . 

(1  +  Y)  |_(Q2/2)  +  g  z\  .  + - *J- 


=  0, 


so,  if  s  is  distance  along  a  streamline,  and  since 


(A. 14) 


J 

(a  +  p) 

- —-11  -  Y  (L  -  c  T)  -  (c  T  +  Y  c  T)  *  0. 

P  ,s  0  p  p  »S 


(A. 15) 

(A. 16) 


If  it  Is  recalled  that 

a  =  I( T)  R'  T,  p  =  p/R  T,  Y  =  Z(T)  R  T/p, 


(A. 17) 


the  last  equation  becomes,  if  p  =  (p  +  a)  and  a  =  (m  /m  ), 

V  a 


R  T 


d  T 
d  p 


♦  [l 

P  -  KT)  L  0 


(c  - 


CP)T 


] 


a  U T) 


[p  -  E(T)]2 


O  Z(T)  cl 
c_  + - E 

+  [l  -  (c  -  c*)T| 

1 

»— 

W 

1 

1  CL 

CL 

_ 1 

L  o  p  J 

a  p  Z'(T) 
£p  -  S(T)]2 


(A. 18) 


where  I'(T)  denotes  the  first  derivative,  and  all  specific  heats  have 

been  taken  as  constant.  In  more  approximate  derivations,  p  =  p, 

(c  +  Y  cv)  *  c  ,  and  c  *  c^  In  the  first  term  of  the  denominator. 

P  P  P  P  P 
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APPENDIX  B.  THE  LOW-LEVEL  INFLOW  LAYER 


The  nature  of  the  incompressible  axisymmetric  flow  under  a  potential  vortex, 
with  negligible  role  for  Coriolis  acceleration,  is  crucial  for  establishing  the 
basis  of  an  inviscid  treatment  of  the  separated  layer.  While  the  general  charac¬ 
ter  of  the  boundary  layer  already  has  been  presented  (Carrier  1971;  McWilliams 
1971;  Burggraf  et  al  1971),  still  a  simpler  and  more  accessible  argument  for 
the  laminar  (constant-viscosity)  case  seems  worth  development  in  this  appendix, 
which  then  goes  on  to  suggest  that  the  character  of  the  laminar  layer  remains 
intact  if  a  gradient-diffusion  model  is  adopted  to  describe  a  turbulent  inflow 
layer.  The  gradient-diffusion  model  of  transfer  in  a  boundary  layer  at  high 
Reynolds  number  is  regarded  here  simply  as  a  convenient  construction  of  empiri¬ 
cal  utility,  without  necessarily  any  basic  validity. 

Previous  work  states  that  the  boundary  layer  grows  monotonical ly,  but 
relatively  slowly,  in  thickness  as  one  proceeds  radially  inward  to  the 
radius  of  maximum  winds  rj  from  the  outer  edge  of  the  closed  system  rQ,  where 
rQ  =  0(10  r j )  suffices  for  most  purposes.  (The  use  of  super  asterisks  to 
denote  dimensional  quantities  is  dispensed  with  in  this  appendix.)  Near  rQ, 
friction  plays  a  role  across  the  entire  vertical  extent  of  the  surface  inflow 
layer;  as  one  approaches  the  axis,  the  role  of  friction  in  the  conservation  of 
radial  momentum  is  confined  to  an  ever-thinning  swirl -free  sublayer  (immediately 
contiguous  to  the  wall)  across  which  the  radial  velocity  component  decreases  to 
zero.  In  accordance  with  the  no  slip  boundary  condition  at  the  wall.  The  pro¬ 
file  of  the  azimuthal  velocity  component  is  shaped  significantly  by  friction 
across  the  entire  thickness  of  the  surface  inflow  layer,  at  all  radial  positions. 


A  nonosci liatory  profile  with  distance  from  the  wall  holds  for  both  radial  and 
azimuthal  velocity  components:  the  maximum  value  of  the  radial  velocity  com¬ 
ponent  can  become  as  large  as  the  potential -vortex  swirl  speed  at  the  same 
radial  distance  from  the  axis  of  rotation,  but  at  fixed  altitude  within  the 
inflow  layer  the  circulation  is  smaller,  at  smaller  radial  position.  Although 
the  swirl  within  the  inflow  layer  sometimes  is  credited  with  values  appreciably 
exceeding  the  asymptotic  value  holding  above  the  inflow  layer  at  the  same 
radial  position  (Lewellen  1976),  such  behavior  is  not  anticipated  here. 

B.l  Formulation 

If,  in  cylindrical  polar  coordinates  (r,0,z),  the  corresponding  velocity 
components  are  denoted  (u,v,w),  the  conservation  of  mass,  radial  momentum,  and 
angular  momentum  may  be  expressed  as  follows  (with  subscripts  r,z  denoting  par¬ 
tial  differentiation  and  v  denoting  kinematic  viscosity): 


(ru)r  +  (rw)z  =  0  ; 


Lr»)»  -  (rv)2  -  (n.)»  .  [v(ru)z]^  -  U(ru)f  -  w(ru)z  ; 
[v(rv)J  -  u(rv)f  -  w(rv)z  *  0  . 


The  boundary  conditions  are  (r  given  const.) 


z  -*■  «:  v  -*■  T/r  ,  T  const  ;  u  -*■  0  ; 


z*0:  u  a  v  ■  w  *  0  ; 


r  ■  r  :  u  *  0,  v  specified  (a  r/r ) 
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The  formulation  is  to  be  applied  only  in  the  range  r^  rQ,  where  r-|  is  the 
(finite)  radius  of  maximum  winds.  Also, 


T  a  r-j  V,  const.  , 


(B.1.7) 


where  V  Is  the  peak  swirl  speed  in  the  potential  vortex.  It  may  be  noted  that 
the  cyclostrophlc  balance  In  the  potential  vortex  [p  denotes  density  and  p 
denotes  pressure] 


p 


(B.1.8) 


has  been  used  to  Identify  the  radial  pressure  gradient  in  (B.1.2). 

B.2  Constant-viscosity  Solution 

If  v  is  constant,  solution  to  (B.1.1)  -  (B.1.7)  is  sought  In  the  forms 
[the  function  form  g(n)  clearly  is  unrelated  to  the  gravitational  acceleration] 

ru  =  -  TB(r)g'(n)  ,  rv  =  Tf(n)  ,  w  *  W(n)  ,  (B.2.1) 

where 


n  = 


r_ 
i  v 


[s(r)]*s 


(B.2. 2) 


Here,  s(r)  and  B(r)  are  functions  to  be  Identified  in  the  course  of  solution. 
One  does  not  expect  (B.2.1)  -  (B.2.2)  to  give  an  exact  solution,  but  to  capture 
the  essence  of  the  functional  behavior  in  a  useful  form. 
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7^- 


Substitution  of  (B.2.1)  and  (B.2.2)  in  (B.1.1)  gives 
wn  .  (vr)®1 

-  (vr)15  g)']i 

“  .  (vr)*  [fi  +  -6C]  g  -  -§C  ng7]  . 
LI  r  ZrsH)  2r$h  J 

where  by  (B.1 .5) 

g(o)  =  o. 


(B.2.3) 


(B.2.4) 


Substitution  of  (B.2.1)  -  (B.2.2)  in  (B.1.2)  and  (B.1.3)  gives,  respec¬ 
tively, 

7T  (1  -  f2  -  B2  g'2)  =  -  Br2  j^l-  -i- +  1^-99'']  -  -£-3B'g'2  ;  (B.2.5) 


*  [4  ■  K  •  *3  *•]  •  •  - 


(B.2.6) 


The  terms  In  (B.2.5)  derive,  from  left  to  right,  from  similarly  positioned  terms 
in  (B.1.2).  For  selfconsistency. 


~  +  7F  -  1  * 

say;  also.  It  Is  necessary,  but  not  sufficient,  to  take 

g'2  *  1  -  f2  , 

Ml  ,  .  1-B2 
r  r* 


From  (B.1.6),  B(rQ)  ■  0;  hence,  from  (B.2.9), 


(B.2.7) 


(B  .2.8) 


(B.2.9) 


r • . 


wen 


l  iJLt  |j 


r  2  -  r2 


Equation  (B.2.7)  may  be  written 


($2s)  =  -  2r$  =»  s  =  —  rQ  (rfl2 


r2)h  , 


(B.2.10) 


(B.2.11 ) 


by  use  of  (B.2.10) 


The  term  In  square  brackets  In  (B.2.5)  Is  0(Br2/s),  while  the  other  terms 

are  0(r2/r2);  substitution  of  (B.2.10)  and  (B.2.11)  Into  these  expressions 

shows  that  the  ratio  of  the  bracketed  to  unbracketed  terms  is  0[(3r2)/(2ro2)], 

so  that  for  r<<  r  ,  the  bracketed  term  is  negligible.  The  neglect  of  the 

bracketed  term  gives  an  Inaccurate  result  where  the  motion  is  slow  [i.e.,  where 

v  «  0(V/10),  u  <  0( V/l 0)].  and  an  accurate  result  where  the  motion  is  fast. 
u  i/ 

From  (B.2.11),  s  n#  (rQ  -  r)  /4  as  r  -*■  rQ;  this  Is  consistent  with  the  more  metic¬ 
ulous  analysis  of  Burggraf  et  al  (1971). 


The  radial  momentum  equation  for  rQ2>>  r2,  i.e.,  for  $  1,  becomes  a 

statement  that  the  left-hand  side  of  (B.1.2)  vanishes;  i.e.,  V2  *  (u2  +  v2), 
such  that  the  radial  component  of  momentum  Is  inviscidly  controlled.  The  two- 
point  nonlinear  third-order  boundary-value  problem,  from  (B.2.5)  -  (B.2.8), 


becomes 


f"  +  gf'  *  0  ; 


(B.2.12) 


1  -  f2 


(B.2.13) 


f(0)  *  g(0)  =  0  ,  f(®)  1 


(B.2.14) 


The  boundary  conditions  follow  from  (B.1.4),  (B.1.5),  (B.1.7),  (B.2.1),  and 
(B.2.2).  The  problem  Is  solved  numerically;  one  finds 


f'(0)  *  0.7456,  g(»)  *  1.941 


(B.2.15) 


•  «.  •  *  •  *,  ■ ,  *  _  '«  * .  .  •  .  1  .  ■  . 


Results  are  graphed  in  figure  30,  where  it  is  evident  that  f  monotonical ly 
increases  (so  g'  monotonical  ly  decreses). 


Of  course,  the  reduction  to  third  order  results  from  the  discarding  of  the 
frictional  term  in  the  conservation  of  radial  momentum  (B.2.5),  which  was  taken 
in  the  inviscid  form  (B .2.8) .  Enforcement  at  all  r  of  u  =  0  at  z  =  0,  which 
implies  g'(0)  »  0,  requires  restoration  of  the  frictional  term  in  a  thin,  effec¬ 
tively  swirl-free  sublayer  contiguous  to  the  boundary. 

A  brief  sketch  of  the  sublayer  across  which  u  falls  to  zero  is  now  pre¬ 


sented.  If 


ru  -  -  r3(r)h'(a)  ,  0  =  -j  ,6  =  1-  (r/rQ)2  , 


then,  from  continuity  (B.1.1), 


(rw)0  “  (t-)5*  m  re,h/  "  re  ~m“  ah^  * 


(B.2.16) 


(B.2.17) 


rw  =  (vr)1*  rnjg'h  -  0  -jjj-  ah'  -  hjj  . 


(B.2.18) 


Substitution  of  (B.2.16)  and  (B.2.18)  in  (B.1.2)  gives 


_ !_.  fjilh'2  .  f  1  lf2 

r*  [raJn  (T^J 


where  the  last  term  involving  f2  Is  dropped  as  uniformly  negligible  over  the 
sublayer.  The  function  m  is  assigned  by  demanding 


(B.2.19) 


(B.2.20) 


(B.2.21) 


If  attention  Is  confined  to  (r/rQ)<<  1,  so  that  8  -*•  1,  8'  =  o(l),  then 

m  =  r  , 
and 

-  h"'  +  (hh')'  =  1,  (B.2.22) 

where  all  neglected  terms  are  0(r2/ro2)  or  smaller.  The  boundary  conditions 
are 

h(0)  =  0,  h'(0)  =  0,  h'H  -*•  1  ;  (B .2.23) 

these  enforce  (B.1.5)  so  that  u(x,0)  =  w(x,0)  =  0,  and  permit  matching  of  fric¬ 
tional-sublayer  and  inviscid-outer-layer  solutions.  One  arrives  at  a  Riccati 
equation: 

-  hf  +  — !j—  =  ~y~  -  h"(0)o  ,  h(o  a>)  -*■  a  ,  (B.2.24) 

where  the  boundary  conditions  (B.2.23)  have  been  used  to  assign  constants  of 
Integration.  Figure  31  shows  that 

h(o)  *  a  +  exp(-o)  -  1  h'(o)  *  1  -  exp(-a)  ,  (B.2.25) 


approximates  the  result  of  numerical  integration  satisfactorily  for  most  pur 
poses.  Hence,  a  uniformly  valid  expression  for  r2 < <  rQ2  is  given  by 


ru 


(B.2.26) 


where  h'(a)  is  given  by  (B.2.24),  or  more  facilely  by  (B.2.25),  and  h'(n)  is 
given  by  (B.2.12)  -  (B.2.14). 


B.3  The  Turbulent  Boundary  Layer  under  a  Rapidly  Swirling  Vortex 

Solution  to  (B.1.1)  -  (B.1.6)  is  again  sought  in  the  forms  (B.2.1),  but 


instead  of  (B.2.2), 


n  =  sTrl  * 


and  B(r),  s(r)  are  now  sought  in  forms  appropriate  for  a  viscosity  taken  in 
the  form 

v  *  —  H(n)  ,  (B 

where  H(n)  is  to  be  specified  and  (B.3. 2)  encompasses  eddy-viscosity  forms  to 
be  examined  here.  It  is  emphasized  at  the  outset  that  selfsimilarity  is  not 
entirely  adequate;  nonselfsimilar  contributions  are  not  uniformly  negligible 
and  warrant  retention  in  improvements  on  the  preliminary  development  (given 
here  to  suggest  that  standard  turbulent  models  for  boundary  layers  yield  results 
with  the  general  characteristics  of  results  already  evolved  for  v  constant). 


Substitution  of  (B.2.1)  into  (B.1.1)  under  (B.3.1)  yields 


-  r  flfii  ♦  Zf]  g  -  al  „9 


Substitution  of  (B.2.1)  into  (B.1.3)  under  (B.3.1)  and  (B.3. 2)  yields 


~r  {[H(n)f ]  '  -  (es)'gf'}  =  0  ; 


consistency  leads  to  the  choice 


(8s)'  ■  -  1  *»  6s  =  r0  -  r  , 


under  (B.1.6).  Hence, 


(Hf')'  +  gf'  -  0  . 


A 


Substitution  of  (B. 2.1 )  and  (B.3.1)  into  (B.1.2)  yields 

[(«  3")'  ♦  gg’]  +  ^  (g ' ) 2  -  |l  -  f!-  e2(g')£]  .  (B.3.7) 

It  is  anticipated  that  the  boundary  layer  associated  with  the  azimuthal 
velocity  component  v  thickens  faster  in  the  turbulent  case  than  in  the  laminar 
case.  The  radial -velocity-component  sublayer,  in  which  friction  enters  sig¬ 
nificantly,  is  anticipated  to  occupy  even  less  of  the  swirl -boundary-layer  thick¬ 
ness  than  in  the  laminar  case.  Thus,  in  most  of  the  inflow  layer,  (B.3.7)  is 
taken  as 

(g')!  *  -  TT  f  -  6!(g’)!]  .  (6.3.8) 


(B.3.9) 

(B.3.10) 

(B.3.11 ) 


It  is  anticipated  that  ultimately  the  definition  of  s  must  be  modified  from 
that  given  In  (B.3.11)  to  account  for  nonnegl Iglble  contributions  from  nonself- 
similar  terms  in  the  equations  for  conservation  of  radial  and  azimuthal  momentum. 


Hence  the  boundary- value  problem  becomes  (B.3.6)  and  (B.3.9),  with  the 
boundary  conditions 


f(0)  *  g(0)  =  0  ;  f(n  -*■  »)  •*  1  . 


(B.3.12) 


It  is  adopted  that  the  total  viscosity  v  -  (vQ  +  vfi),  where  vQ  is  the  molec¬ 
ular  viscosity  and  vg  is  the  eddy  viscosity  and  vg  is  taken  in  the  standard 
form  (Fendell  1972;  Bush  and  Fendell  1972;  White  1974)  (figure  32) 


kl  z2  <uz2  +  vz2^  * 


in  the  wall  layer 


in  the  defect  layer; 


(B.3.13) 


k-|  is  the  square  of  the  von  Karman  const.,  or  (0.41  )2 ;  k2  is  the  Clauser  const., 
0.016;  and  6*  is  the  angular-momentum-defect  thickness.  A  van  Driest-type  wall- 
layer  correction  factor  and  an  intermittency-type  defect-layer  correction  have 
been  omitted  as  nonessential.  In  similarity  variables,  in  view  of  (B.3.9),  (B.3.13) 
becomes 


“tH  k.  ti2  - - — r  Cl  +  f2(02  -  l)]15!  wall  layer 

•  r  >  1  (1  -  f2P 


Ts]  , 

T  k  2  °  • 


(B.3.14) 


defect  layer. 


In  the  defect  layer,  if  o  denotes  a  constant  (to  be  determined). 


6*  -  1  -  -W-  dz  ■  s  (1  -  f)  dn  =  so  . 


(B.3.15) 


Hence,  from  (B.3.2),  for  r<<  rQ  so  6  1 , 


v  *  H(n)  ,  H(n)  =  e  +  k1 


n2f'(n) 


[  [1  -  f2(n)]3s  J 


(B.3.16) 


where  by  definition  subscript  asterisk  here  defines  a  functional  behavior: 


|  F(n)  |  3  F(n)  for  o  3  n  ^  h]  ;  j  F(n)  j  3  F(rV|)  for  n-j  '<  n  '<  00  .  (B.3.17) 


Since  s  rQ  for  r  <  <  rQ  , 


(B.3.18) 


For  consistency  with  similarity,  r  in  (B.3.18)  has  been  taken  uniformly  con¬ 
stant  at  Its  smallest  value  within  the  range  of  interest.  The  eddy  viscosity 
is  made  continuous  by  enforcing  the  following  condition  at  the  Interface 
between  the  viscous  sublayer  (in  which  the  mixing-length  form  holds)  and  the 
defect  layer  (in  which  the  product  of  local  freestream  speed  and  an  integral 
scale  form  a  diffusion  coefficient): 


fr(m) 

k,n?  - n  *  <7  ,  (B .3.19) 

1  ’  [l  -  rtn,)]*1  2 

a  relation  that  defines  rij . 

In  summary,  the  solution  for  the  functions  of  (B.2.1)  and  (B.3.1),  with 
6  identified  in  (B.3.10)  and  s  in  (B.3.11)  and  W  in  (B.3.3),  can  be  obtained 
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from  (B.3.6)  and  (B.3.9)  subject  to  (B.3.12),  where  the  given  parameter  e(  «  1) 
is  defined  by  (B.3.18),  the  derived  parameter  a  is  defined  by  (B.3.15),  and  the 
derived  parameter  is  given  by  (B.3.19),  with  k^=(0.41)2  and  k£  =  0.016. 

Numerical  solution  entails  a  two -degree-of -freedom  iteration  in  that,  for 
example,  one  must  anticipate  both  f'(0)  and  n-j ,  and  confirm  the  validity  of  the 
guesses  by  a  compatible  with  (B.3.15)  and  f(n  ®)  compatible  with  (B.3.12),  in 
a  shooting-type  approach.  Facile  reduction  to  a  one-degree-of-freedom  iteration 
is  not  readily  executed.  Nevertheless,  solutions  for  various  values  of  e  are 
straightforward,  and  are  presented  in  figures  33  and  34,  with  the  values  required 
for  "shooting"  tabulated  in  table  2. 
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FIGURE  1.  Schematic  of  the  postulated  four-part  model  of  the  structure 

mature  two-cell  vortex,  of  radial  extent  r*  and  axial  extent  zi  (the 

0  » 

height  of  the  tropopause). 
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FIGURE  2.  Schematic  of  the  location  of  an  inviscid  eyewall  without  structure 
demarcating  the  interface  between  an  isobaric  nonswirling  eye  at  pressure 
P|ye  an<*  a  Potent^  vortex.  The  sheet  representing  the  eyewall  has  thickness 
h*  0;  in  the  turnaround  region,  its  displacement  from  the  axis,  as  a  function 
of  height  above  the  ground  plane,  is  denoted  R*(z*). 
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FIGURE  3.  Phase-plane  properties  of  equation  (3.1.11),  where  P  =  (dR/dz), 
with  R  the  dependent  variable  (R  >  0)  and  z  the  independent  variable. 
Isoclines  of  zero  and  infinite  slope  are  noted.  At  R  *  1 ,  the  distance 
from  the  axis  of  symmetry  at  which  the  surface  Inflow  layer  separates, 
a  finite  positive  slope  is  adopted.  For  a  >  0,  R*  >  1;  for  a  =  0,  R,  !  1 
The  sketched  trajectory  (solution  curve)  Is  a  limit  cycle  (closed  curve 
Indicative  of  periodic  behavior).  The  periodicity  is  not  of  physical 
Interest. 


solution  to  the  boundary- value  problem  posed  by  equations  (3.1.7) 
,  =  2;  the  solution  may  be  completed  by  symmetry  considerations  to 


denotes  the  value  of  (1  -  R-2);  the  "swirl”  denotes  the  value  of  aR”3(l  +  R'2)-^;  the  quantity 
"radlal/axlal"  denotes  the  value  of  "swirl"  minus  "pressure", i.e. ,  the  value  of  the  second 
term  on  the  right-hand  side  of  (3.1.7)  before  solution  for  the  second  derivative. 
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FIGURE  6b.  Pressure- temperature  relationships. 


FIGURE  7.  The  radial  displacement  of  the  outer  edge  of  the  eyewall  R. 
i.e.,  the  "ambient  edge"  contiguous  to  the  potential 
vortex,  as  a  function  of  height  above  ground  level  z,  for 
a  tornado- like  case. 


FIGURE  8.  For  the  case  of  figure  7,  the  separated- layer  thickness  h, 
the  transversely  averaged  swirl  speed  v,  the  transversely 
averaged  radial-axial-flow  speed  q,  and  the  pressure  at 
the  ambient  edge  of  the  separated  layer  pag,  as  functions  of 
of  altitude  above  ground  level  z.  The  pressure  in  the  eye 
1s  he1d  at  unity. 
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FIGURE  10.  For  the  case  of  figure  9,  the  separated- layer  thickness  h, 
the  transversely  averaged  swirl  v,  and  the  transversely 
averaged  radial-axial-flow  speed  q,  as  functions  of  alti¬ 
tude  above  sea  level  z. 
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FIGURE  11.  For  the  case  of  figure  9,  the  pressure  at  the  ambient  edge 
of  the  separated  layer  pfle  and  the  pressure  in  the  eye  Peye» 
as  functions  of  altitude  above  sea  level  z. 


Highly  convectlvely  unstable  ambient  (measured  in  the  vicinity  of 
a  tornado)  at  St.  Cloud,  MN,  on  July  22,  1967  at  1715Z. 
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FIGURE  13.  Pressure* temperature  relationships  for  the  ambient  of 

figure  12,  for  ground-level  air  raised  on  a  moist  adiabat 
to  its  level  of  neutral  stability  ( 1 . e. ,  to  the  tropo- 
pause  as  defined  here),  and  for  tropopause-level  air  com¬ 
pressed  dry-adiabatically.  Heights  are  associated  by 
means  of  hydrostatics,  with  z*=0  being  ground  level.  The 
fractional  heights  displayed  are  the  result  of  conversion 
from  engineering  units  to  SI  convention. 
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FIGURE  14.  For  the  case  of  figure  12  with  a  fully  Inserted  eye,  the 
adopted  radial  displacement  R  of  the  outer  eyewall  ( 1 . e . , 
the  eyewall  edge  contiguous  to  the  potential  vortex)  as  a 
function  of  altitude  z,  where  the  completion-of-turnaround 
height  Is  taken  as  z  *  0.  The  Inner  eyewall  position  Is  the 
site  of  the  eyewall  edge  contiguous  to  the  eye;  the  Inner 
edge  Is  located  by  subtracting  the  eyewall  thickness  h(z)  In 
a  direction  perpendicular  to  the  local  tangent  to  R(z). 
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FIGURE  15.  For  the  case  of  figure  12  with  a  fully  inserted  eye,  the 
eyewall  thickness  h,  and  the  transversely  averaged  eyewall 
radial -axial -flow  speed  q,  as  functions  of  altitude  z,  for 
the  R(z)  adopted  in  figure  14. 
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FIGURE  17.  For  the  case  of  figure  12  with  a  fully  inserted  eye,  the 
angular  momentum  r  as  a  function  of  altitude  z,  for  the 
R(z)  adopted  in  figure  14.  The  "jitter"  is  attributed  to 
simplistic  numerical  procedure. 
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FIGURE  18.  For  the  case  of  figure  12  with  a  fully  inserted  eye,  the 
transversely  averaged  eyewall  swirl  speed  v,  and  the  swirl 
speed  of  the  potential  vortex  at  the  eyewall  edge  (r/R),  as 
functions  of  altitude  z.  Maintaining  (r/R)>v  until  high 
In  the  troposphere  constrains  the  choice  of  R(z)  adopted 
(figure  14). 
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FIGURE  19.  For  the  case  of  figure  12  with  a  90%  inserted  eye,  the 
adopted  radial  displacement  R  of  the  outer  eyewall  as  a 
function  of  altitude  z. 
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FIGURE  21.  For  the  case  of  figure  12  with  a  90%  Inserted  eye,  the 
angular  momentum  r  as  a  function  of  altitude  z,  for  the 


R(z)  adopted  In  figure  19. 


FIGURE  22.  For  the  case  of  figure  12  with  a  90%  inserted  eye,  the  eyewall 
swirl  speed  v,  and  the  swirl  speed  of  the  potential  vortex 
at  the  eyewall  edge  (r/R),  as  functions  of  altitude  z.  Satis¬ 
faction  of  the  stability  constraint  (r/R) >  v  through  most  of 
the  troposphere  furnishes  a  strong  constraint  on  the  iteratively 
arrived-at  selection  of  R(z)  presented  in  figure  19. 
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FIGURE  23.  Moderately  convectlvely  unstable  ambient  (measured  In  the 
vicinity  of  a  tornado)  at  Jackson,  MS,  on  April  17,  1978  at 
2300Z. 
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FIGURE  24.  Pressure- temperature  relationships  for  the  ambient  of  figure  23, 
for  ground-level  air  raised  on  a  moist  adlabat  to  Its  level  of 
neutral  stability  (l.e.,  to  the  tropopause  as  defined  here),  and 
for  tropopause- level  air  compressed  dry-adlabatlcally.  Heights 
are  associated  by  means  of  hydrostatics,  with  z**0  being  ground 
level.  The  fractional  heights  displayed  are  the  result  of  con¬ 
version  from  engineering  units  to  SI  convention. 
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FIGURE  25. 


For  the  case  of  figure  23  with  a  fully  Inserted  eye,  the 
adopted  radial  displacement  R  of  the  outer  eyewall  (l.e., 
the  eyewall  edge  contiguous  to  the  potential  vortex)  as  a 
function  of  altitude  z,  where  the  completlon-of- turnaround 
height  Is  taken  as  z*0.  The  Inner  eyewall  position  Is  the 
site  of  the  eyewall  edge  contiguous  to  the  eye;  the  Inner 
edge  Is  located  by  the  eyewall  thickness  In  a  direction 
perpendicular  to  the  local  tangent  to  R(z). 
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FIGURE  26.  For  the  case  of  figure  23  with  a  fully  Inserted  eye,  the 
eyewall  thickness  h,  and  the  transversely  averaged  eyewall 
radlal-axlal-flow  speed  q,  as  functions  of  altitude,  for 
the  R(z)  adopted  In  figure  25. 


K 


FIGURE  27.  For  the  case  of  figure  23  with  a  fully  Inserted  eye,  the 


pressure  In  the  eye  p  and  the  pressure  at  the  ambient 

eye 


edge  of  the  eyewall  pafi  as  functions  of  altitude  z,  for 


the  R(z)  adopted  In  figure  25. 
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FIGURE  28.  For  the  case  of  figure  23  with  a  fully  Inserted  eye,  the 
angular  momentum  r  as  a  function  of  altitude  z,  for  the 
R(z)  adopted  In  figure  25.  The  "jitter"  Is  attributed  to 
simplistic  numerical  procedure. 
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FIGURE  29.  For  the  case  of  figure  23  with  a  fully  inserted  eye,  the 
transversely  averaged  eyewall  swirl  speed  v,  and  the  swirl 
speed  of  the  potential  vortex  at  the  eyewall  edge,  as  func¬ 
tions  of  altitude  z.  Maintaining  (r/R)  >  v  until  high  in 
the  troposphere  constrains  the  choice  of  R(z)  adopted 
(figure  25). 
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'eloclty  component.  The  similarity  Independent  variable  n  is  large  for  large  distance 
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component  v,  and  axial  velocity  component  w,  at  a  radial  distance  r  from  the  center  of  the  vortex 
at  which  the  impressed  swirl  V  is  large. 


0.40 


IRE  33.  Numerical  results,  for  Inverse  Reynolds  number  i 
turbulent  boundary  layer  under  a  rapidly  swirling  vortex 
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0.550 

1.709 

1.772 

0.2 

0.4 

O.S 

10.0 

0.S94 

0.S2S 

1.030 

1.032 

0.2 

0.4 

O.S 

100.0 

0.570 

0.509 

1.077 

1.007 

TABLE  1.  Numerical  results  (denoted  "calc")  for  extrema  of  the  turnaround 
from  integration  of  the  Initial-value  problem  (3.1.7)-(3.1.9).  Results  for 
are  ^rom  (3.1.16).  Also,  calculated  results  for  R*  for  RQ'  =100  are 
within  2%  of  results  for  RQf  +  *>,  for  a, 8  studied. 


*4  " 


e 


f'(0) 


nl 

10"*  1.35*10"2  6.18*10-*  1. 17868* 102 

10- 5  6.90*10-*  3.85*10"*  6.24682*102 

10"7  2.88*10-*  2.395*10"*  2.56536*104 

10"*  1.637*10-*  1.805*10-*  1. 36583*10® 

10"lx  1.034*10"*  1.436*10"*  9.72465*107 

TABLE  2.  Properties  of  the  approximate  selfsimilar  solution  for  the 
turbulent  boundary  layer  under  the  high-speed  portion  of  a  potential 
vortex. 


